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SUMMARY 


An  analytical  and  experimental  investigation  was  conducted  to  acquire 
systematic  model  rotor  performance  and  wake  geometry  data  and  to  evaluate 
the  accuracy  of  various  analytical  methods  in  predicting  the  effects  on 
performance  of  changes  in  helicopter  rotor  design  and  operating  parameters. 
Both  classical  hover  performance  analyses  and  analytical  methods  recently 
developed  at  the  United  Aircraft  Research  Laboratories  were  evaluated. 

Of  primary  concern  in  the  study  was  the  assessment  of  assumptions  in  the 
analyses  regarding  the  geometry  of  the  rotor  wake.  It  was  found  that 
analyses  based  on  a  contracted  wake  geometry  generally  provided  signifi¬ 
cantly  improved  predictions  of  performance  for  those  rotor  operating 
conditions  where  the  more  classical  uncontracted  wake  analyses  exhibited 
major  shortcomings.  Attempts  to  develop  a  theoretical  method  for 
predicting  contracted  wake  geometries  were  only  partially  successful 
although  the  method  yielded  good  qualitative  results.  Of  particular 
interest  was  the  prediction  by  the  analysis  of  an  instability  of  the  tip 
vortex  helix  at  moderate  distances  from  the  rotor  which  appeared  to  be 
substantiated  by  available  experimental  results.  In  lieu  of  an  accurate 
theoretical  wake  method,  it  was  recommended  that  the  experimental  wake 
data  measured  in  this  study  provide  the  contracted  wake  geometries  needed 
in  performance  calculations.  Analysis  of  the  wake  data  for  that  portion  of 
the  wake  which  was  stable  (i.e.,  near  the  rotor)  indicated  that  the  data 
could  be  expressed  in  relatively  simple  generalized  equations  which 
facilitate  the  rapid  estimation  of  contracted  wake  geometries  for  a  wide 
range  of  rotor  designs  and  operating  conditions.  Finally,  it  was 
demonstrated  that  rotor  performance  is  sensitive  to  small  changes  in  the 
position  of  the  tip  vortex  relative  to  the  following  blade,  and  it  was 
recommended  that  additional  full  scale  correlation  studies  be  made  to 
provide  further  information  on  the  adequacy  of  the  generalized  wake 
geometry  charts  provided  herein. 
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INTRODUCTION 


The  need  for  attaining  peak  lift  system  performance  is  greater  with 
rotary-wing  VTOL  aircraft  than  with  conventional  aircraft.  This  results 
directly  from  the  generally  lower  payload  to  gross  weight  ratio  of  such 
aircraft  which,  in  turn,  increases  the  payload  penalty  associated  with 
any  unexpected  deficiencies  in  performance  that  might  arise  as  a  result  of 
shortcomings  in  the  design  analyses  employed.  For  example,  since  the  pay- 
load  is  typically  2%  of  the  gross  weight,  a  performance  deficiency  of  1% 
in  lift  capability  can  result  in  a  reduction  in  payload. 

As  described  in  Reference  1  .mmonly  used  theoretical  methods 
bec-ane  inaccurate  as  number  of  blades,  blade  solidity,  blade  loading,  and 
tip  Mach  number  are  increased.  The  discrepancies  noted  appear  to  stem 
from  simplifying  assumptions  made  in  the  analyses  regarding  the  geometric 
characteristics  of  the  rotor  wake.  In  Reference  1,  a  method  for 
considering  the  effects  of  wake  contraction  on  hover  performance  was 
introduced.  This  computerized  method,  developed  at  the  United  Aircraft 
Research  Laboratories  (UARL)  and  termed  the  UARL  Prescribed  Wake  Hover 
Performance  Method,  requires  a  prior  knowledge  of  the  wake  geometry. 
However,  at  the  time  Reference  1  was  written  (1967),  available  wake 
geometry  data  were  extremely  limited.  Due  to  the  expense  involved, 
systematic  wake  geometry  data  on  full-scale  rotors  were  almost  nonexistent. 
Available  model  results,  on  the  other  hand,  were  limited  to  rotors  having 
three  blades  or  less  and  operating  at  low  tip  Mach  numbers.  Thus,  two 
methods  of  approach  were  initiated  under  this  investigation  to  obtain  the 
required  wake  geometry  information.  In  the  first  an  experimental 
investigation,  using  model  rotors,  was  conducted  in  which  a  systematic, 
self-consistent  set  of  data  on  rotor  performance  and  associated  wake 
geometry  characteristics  was  obtained  for  a  wide  range  of  blade  designs 
and  operating  conditions.  In  the  second,  an  available  analytical  method 
for  predicting  rotor  wake  geometry  in  forward  flight,  describe  1  in 
Reference  2,  was  extended  to  the  hover  condition.  Briefly,  the  method 
developed  involves  the  establishment  of  an  initial  wake  model  comprised 
of  finite  vortex  elements  and  the  repeated  application  of  the  Biot-Savart 
law  to  compute  the  velocity  induced  by  each  vortex  element  at  the  end 
points  of  all  other  vortex  elements  in  the  wake.  These  velocities  are 
then  integrated  over  a  small  increment  of  time  to  determine  the  new 
positions  of  the  wake  elements,  and  the  entire  process  is  repeated  until 
a  converged  wake  geometry  is  obtained. 
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The  incorporation  of  the  experimental  and  analytical  wake  geometry 
in  the  Prescribed  Wake  Method  results  in  two  analyses  (the  described 
Experimental  Wake  Analysis  and  the  described  Theoretical  Wake  Analysis) 
for  computing  hover  performance.  The  availability  of  model  rotor  data 
permits  the  evaluation  of  these  analyses  by  (l)  providing  experimental 
wake  data  both  for  input  to  the  described  Experimental  Wake  Analysis  and 
for  comparison  with  predicted  wake  geometry  results  of  the  Prescribed 
Theoretical  Wake  Method,  and  (2)  providing  consistent  experimental  perfor¬ 
mance  data  for  comparison  with  predicted  performance  results.  Thus,  the 
principal  objectives  of  this  investigation  were  to: 

(a)  Provide  experimental  information  on  the  performance  and  wake 
geometry  characteristics  of  hovering  model  rotors  as 
influenced  by  number  of  blades,  blade  twist,  blade  aspect 
ratio,  rotor  tip  speed,  and  blade  collective  pitch  setting 

(b)  Modify  an  existing  forward-flight  distorted  wake  program  to 
permit  the  prediction  of  the  wake  geometry  characteristics 
in  hover 

(c)  Evaluate  the  accuracy  of  various  hover  performance  theories 
having  differing  rotor  wake  geometry  assumptions 

Included  in  this  report  are:  (1)  a  description  of  the  model  rotor 
experimental  program,  (2)  a  discussion  of  the  experimental  rotor 
performance  and  wake  geometry  results,  (3)  comparisons  of  the  experi¬ 
mental  wake  geometry  results  with  other  experimental  sources, 

(4)  descriptions  of  the  theoretical  methods  for  predicting  wake  geometry 
and  hover  performance,  (5)  a  discussion  of  the  results  of  the  evaluation 
of  the  wake  geometry  analysis,  and  (6)  a  discussion  of  the  results  of  the 
evaluation  of  the  theoretical  methods  for  predicting  hover  performance. 
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MOEEL  ROTOR  HOVER  TEST 


TEST  EQUIPMENT 
Model  Test  Facility 

The  test  program  vu  conducted  on  a  model  rotor  hover  test  facility 
located  at  the  United  Aircraft  Research  Laboratories  (UARL).  The  test 
facility,  shown  in  Figure  1,  consists  of  a  large  enclosed  area  approxi¬ 
mately  45  by  55  fee t  with  a  ceiling  height  of  40  feet.  The  facility  is 
equipped  with  a  rotor  test  rig,  flow  visualisation  equipment,  and  a 
movable  ground  plane  which  was  considered  to  be  in  an  out-of-ground-effect 
position  when  lowered  to  3.5  rotor  radii  below  the  rotor.  A  whirl-stand 
model  is  also  available  for  the  simulation  of  conditions  on  the  Sikorsky 
Aircraft  full-scale  whirl  stand.  This  facility  is  the  same  one  used  for 
the  model  tests  reported  in  Reference  1. 

A  photograph  of  the  model  rotor  test  rig  is  shown  in  Figure  2.  A 
40-horsepower,  variable-speed  electric  motor  was  used  as  a  power  source. 
The  rotor  was  driven  through  a  3:1  speed  reduction  system  to  allow  opera¬ 
tion  at  a  tip  speed  of  700  ft/sec  at  mart  nun  available  power.  Average 
rotor  thrust  and  torque  measurements  were  made  by  means  of  strain-gaged 
load  cells  mounted  above  the  rotor  an  a  support  frame.  The  motor-balance 
assembly  is  shown  schematically  in  Figure  3.  Additional  instrumentation, 
used  to  monitor  operation,  included  a  solid-state  counter  for  measuring 
rpn,  vibration  maters,  and  a  model  power  control  console. 

Flow  visualisation  equipment  included  the  following: 

1.  Ammonium  sulphite  “smoke". 

2.  Variable-position  smoke  raJtet;. 

3.  Two  high- intensity,  short -duration  light  sources  (ederoflash 
units)  for  stop-action  still  photographs. 

4.  High-intensity  lights  for  high  frame-speed  movies. 

3,  An  electronic  time-delay  control  for  cameras  and  micro  flesh 
unit*  to  permit  systematic  pbotogrephing  of  the  cyclic  time 
history  of  the  rotor  wake. 

6.  Polaroid,  70  mm,  and  Fas  tax  movie  cameras. 
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Model  Rotors 


The  model  rotor  system  consisted  of  a  multibladed  rotor  hub  and 
specially  designed  model  blades.  The  rotor  hub,  shewn  in  Figure  4,  was 
designed  to  accommodate  any  number  of  blades  up  to  and  including  eight. 
Flapping  hinges  were  provided  (but  no  lag  hinges),  and  blade  collective 
pitch  was  varied  manually. 

Four  sets  of  model  rotor  blades  were  used  to  conduct  the  test 
program.  The  model  blade  design  consisted  of  an  aluminum  spar  and  a 
balsa  trailing -edge  section,  as  shown  in  Figure  5.  The  blades  were 
untapered  and  designed  such  that  the  elastic  axis,  section  center  of 
gravity,  and  center  of  pressure  were  coincident  at  the  quarter-chord 
position  (within  1%  of  the  chord).  The  mass  and  stiffness  properties  of 
the  model  blades  greatly  exceeded  those  of  model  blades  dynamically  scaled 
to  typical  full-scale  blades.  For  example,  the  Lock  number  of  the  blades 
with  an  aspect  ratio  of  18.2  operating  at  a  tip  speed  of  TOO  fps  was 
approximately  3.0  compared  to  a  typical  full-scale  Lock  number  of  10. 
Hence,  model  blade  coning  angles  were  lower  than  full-scale  coning  angles. 
However,  the  use  of  such  rotor  blades  permitted  concentration  on  the 
aerodynamic,  rather  than  aeroelastic,  aspects  of  rotor  hover  performance. 
The  blade  parameters  are  summurited  in  Table  I. 


TABLE  I. 

!***!,  ROTOR 

BLADE  PARAMETERS 

Blade  Parameter 

Design  (1) 

Design  (2) 

Design  (3) 

Design  (6) 

Linear  Twist,  B^(deg) 

0 

-8 

-16 

-8 

Aspect  Ratio,  AR 

10.2 

18.2 

18.2 

13.6 

V 

26.75 

26.75 

26.75 

26.75 

chord,  c  (in.) 

1.47 

1.47 

1.4? 

1.96 

Airfoil  Section  (BACA) 

0012 

0012 

0012 

0012 

Root  Cutout  (|R) 

14.8 

14.8 

14.8 

14.8 

j  Flapping  King*  Offset  ($R)  6.8 

6.8 

6.8 

6.8 

1  Taper 

Bone 

a*-—-. — 

non* 

Bone 

The  second  set  of  blades  listed  in  Table  I  (8j  =  -8  deg,  AR  =  18.2)  was 
constructed  and  performance-tested  by  United  Aircraft  prior  to  this 
investigation.  At  the  time  of  this  investigation,  four  of  these  blades 
were  available  in  their  original  form  and  eight  were  available  in  a 
modified  form.  The  modification  consisted  of  the  adaption  of  a  plastic 
tip  section  to  the  outer  1 2%  of  the  blade.  This  tip  section  was  nominally 
identical  in  shape  to  the  removed  aluminum-balsa  section.  The  plastic  tip 
blades  were  used  in  this  investigation  to  obtain  ground-effect  and  flew 
visualization  results  for  the  six-  and  eight-blade  configurations  of  this 
particular  blade  design.  Other  performance  results  for  this  desigr  are 
included  from  the  Ubited  Aircraft  investigation  mentioned  above.  Hie 
blades  corresponding  to  the  three  other  blade  designs  were  fabricated 
specifically  far  this  investigation. 

Each  rotor  blade  design  was  tested  with  the  following  numbers  of 
blades  and  rotor  solidity  ratios: 


Humber  of  Blades,  b 
2  4  6  8 

Rotor  Solidity  Ratio,  ^(Designs  (l)-(3))  0.035  0.070  0.105  0.140 
Rotor  Solidity  Ratic,  <r (Design  (4))  0.0467  0.0933  0.140  0.1867 


TEST  HtOCBDURES 
Calibration 

ftior  to  testing,  the  thrust  and  torque  instrumentation  was 
calibrated.  The  thrust  load  cell  was  calibrated  on  the  test  rig  by 
hanging  weights  free  the  rotor  shaft,  Jbe  rig  was  calibrated  in  torque 
by  suspending  waists  through  e  known  moment  arm.  The  thrust  and  torque 
calibration  derivativas  wart  determined  directly  in  strain  gage  units 
per  pound  (SOUS /lb)  and  strain  gaga  units  par  foot-pound  (SCUS/ft-lb), 
respectively. 

A  dowel  pin,  mounted  perpendicular  to  an  arm  astending  from  the  root 
of  the  blade  along  the  blade  chord  (see  Figure  4),  was  used  to  manually 
set  collective  pitch  angle  relative  to  a  flat  on  the  blade  retention 
fitting  (attached  to  the  hub).  For  each  blade,  the  distance  from  the  pin 
to  the  flat,  measured  by  means  of  a  depth  micromater,  was  calibre  ted  with 
respect  to  the  collective  pitch  angle  at  the  three-quarter  radius.  Slade 
tracking  wet  checked  by  Observation  of  tha  blada  tips  through  a  transit, 
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with  lighting  supplied  by  a  Strdbotac  triggered  at  a  specified  nuaber  of 
flashes  per  rotor  revolution.  In  this  Banner,  several  blades  were 
observed  at  once,  and  their  tip  positions  were  coopered. 

To  calibrate  the  flow  visualisation  photographs  and  to  ainiaize 
errors  due  to  caaera  angle  and  lens  distortion,  a  planar  grid  indicating 
increaents  of  the  rotor  radius  (0.475  in.)  was  placed  in  the  plane  of 
the  saoke  (reference  plane)  and  photographed  prior  to  the  test.  Photo¬ 
graph s  of  this  grid  systea  (Figure  6)  were  used  in  the  construction  of  e 
grid  teaplate  overlay  for  the  reduction  of  the  flow  visualisation  photo¬ 
graphs  to  radial  and  axial  wake  coordinates.  The  blade  axiautb  position 
was  calibrated  for  each  rpa  by  calculating  the  delay  tine  between  the 
passage  of  a  reference  blade  through  the  reference  plane  and  the  passage 
of  a  single  gear  tooth  aounted  on  the  rotor  sLaft. 

Test  Iterant ters 


Systematic  data  were  obtained  to  aeasure  the  effect  of  the  following 
paraaeters  on  rotor  hover  perfcraance  and  associated  wake  gecaetry 
characteristics . 


Prisary  Test  fturaMter 


Mfinal  Test  Values 


1.  Buber  of  BUdaa,  b 

2.  Blade  linear  twist,  9^ 
$.  Blade  aspect  ratio,  AR 

4.  Rotor  tip  speed, HR 

5.  Collective  pitch, 


2,  4,  6,  B 
0,  -B,  -16  deg 
13.6  end  18.2 
525,  600,  TOO  fpe 
0  to  MX* 


•Determined  by  operating  stall  liaits 


Variationa  in  Urea  of  the  above  primary  test  psremt era  were  equi valent 
to  independent  variations  in  three  related  peraanrters.  That  la,  varia¬ 
tions  In  number  of  bladaa  were  equivalent  to  variations  in  rotor  solidity, 
aa  specified  previously.  Variationa  in  rotor  tip  speed  of  $25,  600,  end 
700  fpe  were  equivalent  to  variations  in  tip  ttsch  number  of  0.46,  0.525, 
and  0.6l,  respectively.  Finally,  variations  in  collective  pitch  were 
used  to  produce  variations  in  rotor  thrust  level.  Roadnal  collective 
pitch  settings  of  0,  6,  8,  10,  and  12  deg  were  used  whenever  possible. 
Additional  values  were  tested  to  provide  aore  extensive  data  for  sqm  test 
configurations ,  particularly  In  the  stall  region  where  the  eexieue 
collective  pitch  was  Halted  at  the  higher  tip  speeds  by  an  apparent  stall 
flutter. 
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Data  regarding  the  effect  of  the  parameters  listed  above  were 
obtained  with  rotors  operating  substantially  out  of  ground  effect 
(height  above  the  ground  plane  equal  to  3 .5  rotor  radii).  In  addition  to 
the  out-of-ground-effect  (OGE)  conditions,  a  number  of  data  points  were 
taken  in  ground  effect  (IGE)  and  with  a  model  whirl  stand  simulating  the 
Sikorsky  Aircraft  full-scale  rotor  whirl  stand  shown  in  Figure  7. 


Die  test  paraaeter  coabinations  (test  conf iguraticns )  that  were 
investigated  are  given  in  Table  II. 


TABUS  II .  TEST  PARAMETER  CCMBMATIOHS 


Linear 

Twist 

(deg) 

Aspect 

Ratio 

Rotor 

Condition 

Rotor 

Height/ 

Radius 

Ho.  of 
Blades 

Ho.  of 
Tip 
Speeds 

Ho.  of 
Collective 

Pitch  Values 

0 

18.2 

OGB 

3.5 

2,**  ,6,8 

3 

5* 

-a 

18.2 

0GB 

3.5 

**  — 

2, **,6,8 

3 

5* 

-16 

18.2 

OGB 

3.5 

2,*»  ,6,8 

3 

5* 

13.6 

OGB 

3.5 

2 ,4, 6, 8 

3 

5* 

-8 

18.2 

IGE 

1.67 

2,«t6,8 

3 

1 

-a 

18.2 

IGB 

0.67, 

1*0,1.33, 

1.67,2.0 

6,8 

3 

2 

*8 

18.2 

Whirl 

Stand 

1.67 

2, 4, 6, 8 

3 

1 

•Minimum  nuaber. 

•vAnfanaaoi  results  obtained  tram  previous  investigation. 

Data  Acquisition 


The  procedure  for  data  acquisition  fundt— ntally  consisted  of  setting 
the  test  configuration  (collective  pitch,  ground  plane  position,  and 
number  of  blades)  and  then  recording  perform nee  and  wake  geo— try  data 
at  the  required  tip  speeds.  The  aabient  temperature  and  pressure  in  the 
enclosed  area  were  monitored  and  recorded  during  the  test.  The  rotor 
thrust  and  torque  data  were  obtained  by  Manually  recording  the  outputs  of 
the  thrust  and  torque  load  cells  on  self-balancing  potentio— tore  in 
strain  gage  units  (S0U5).  Oscillations  of  the  potent! a— ter  readings  es 
high  as  *10  SC  US  (<0.03  lb)  for  thrust  and  *30  SOUS  (*0.2  ft, -lb)  were 
observed  far  so—  test  conditions.  To  obtain  representative  steady-state 
values,  average  readings  were  recorded.  Bach  test  condition  wae  repeated 
twice  within  a  test  run  (a  test  run  consisted  of  data  recorded  between  the 
starting  end  stopping  of  the  rotor  rotation),  end  the  results  were 
averaged.  In  addition,  aost  test  conditio—  were  repeated  in  two 
consecutive  runs,  tony  of  the  test  conditio—  were  also  rep— ted  an 
different  dates  to  check  the  repeatability  of  the  date. 

For  aost  rotor  configuratla— ,  the  — 1—  collective  pitch  at  which 
operation  was  possible  was  limited  by  the  occurrence  of  e  rapid  increase 
in  rotor  noise  es  tip  speed  was  increased  beyond  e  certain  level.  A 
strain  gage  to  measure  torsional  response  wee  placed  at  the  root  of  the 
blade  and  substantial  Increases  in  blade  torsional  response  were  noted 
under  these  conditio—,  and  it  wee  inferred  that  e  stall  flutter  houndery 
wes  being  penetrated.  Collective  pitch  (thrust  level)  at  a  given  tip 
speed  was  then  limited  by  what  will  he  deacribed  as  e  stall  flutter 
boundary  in  the  discussion  of  results. 

To  obtain  flow  visualisation  data,  smoke  was  injected  into  the  flaw 
by  movable  smoke  rekes  located  above  end  to  both  sides  of  the  rotor 
(figure  2).  Aemonium  sulphite  smoke  was  generated  by  nixing  ammonium  gas 
and  sulphur  dioxide  gas  from  separate  peats  ©a  each  smoke  nestle  mounted 
on  the  smoke  rekes.  The  smoke  rekes  ware  positioned  by  ran ote  control  to 
insure  a  clearly  defined  tip  vortex.  The  wake  patterns  ware  recorded  on 
filn  with  remotely  operated  cameras.  Illumination  was  provided  for  still 
photographs  by  two  microfltsh  units  (tine  duration  0.3  microsecond  for 
stop  action)  and  for  high  frame-speed  movies  by  sixteen  350-watt  flood¬ 
lights.  A  tine -delay  system  was  used  to  trigger  the  cameras  and  micro-. 
Pash  units  when  the  rotor  was  at  a  desired  asinuth  position.  The  delay 
system  used  a  ona/rev  signal  from  the  rotor  shaft  as  a  reference,  and  the 
delay  ti—  (manually  adjusted)  was  measured  on  an  electronic  counter. 

For  each  test  condition,  70—  photographs  were  taken  at  preselected 
asinuth  positions  af  a  reference  blade  with  respect  Is  the  plane  of  the 
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smoke.  For  example,  for  a  two-bladed  rotor,  photographs  were  generally 
taken  at  azianth  angles  of  0,  15,  30,  60,  90,  120,  and  150  deg.  For  a 
six-bled ed  rotor,  azianth  angles  were  generally  0,  15,  30,  and  45  deg. 
■carnally,  two  or  three  photographs  were  taken  at  each  aziauth  position. 
IWo  caaeras  were  need  to  Obtain  photographs  of  both  the  complete  rotor- 
wake  systea  and  dose- up  views  of  the  right  half  of  the  rotor-wake  system. 
To  supplement  these  still  photographs,  high  fraae- speed  movies 
(2000  frame  s/sec)  were  taken  at  selected  conditions. 


data  mucnm 


Thrust  and  torque  measurements  (in  SGUB)  were  converged  to  thrust  (lb) 
and  torque  (ft- lb),  thrust  coefficient  (Cy)  and  torque  coefficient  (Cq), 
end  thrust  coefficient-solidity  ratio  (Cy/cr)  mad  torque  coefficient- 
solidity  ratio  (Cq/ff)  values  for  all  test  conditions.  The  air  density  (P) 
used  In  nondiaansionalizing  the  data  was  calculated  for  each  test  condi¬ 
tion,  based  on  the  recorded  temperature  and  pressure  readings. 


Flow  Visualisation  Data 


To  introduce  the  procedures  used  in  reducing  the  flow  visualisation 
data,  a  brief  dlaeusslon  of  the  fundamental  characteristics  of  the  hover¬ 
ing  rotor  woke  sad  the  interpretation  of  thaae  characteristics  from  smoke 
photographs  is  presented.  A  schematic  of  the  weke  from  one  blade, 
reproduced  directly  from  (kray's  Interpretation  in  haference  3,  is  shown  in 
Figure  8.  The  weke  contains  two  primary  eosgone-it*.  The  first,  and  most 
prominent,  is  tke  strong  tip  varies  which  arise*  free  the  rapid  roUlt« 
up  of  the  portion  of  the  vortex  sheet  shad  from  the  tip  region  of  the 
blede.  The  ascend  feature  is  the  vortex  sheet  shed  free  the  inboard 
section  of  the  blade.  This  sheet  doss  not  roll  up  but  generally  remains 
in  the  form  of  distributed  verticlty.  tbe  vertical  or  axial  transport 
velocity  near  the  outer  and  of  the  inboard  vortex  sheet  is  such  greater 
than  that  of  the  tip  vortex.  The  vertical  velocity  of  the  inboard  sheet 
also  increases  with  radial  petition,  resulting  in  a  substantially  linear 
cross  section  of  the  inboard  sheet  at  any  specific  azimuth  position  as 
shown  in  figures  9  and  10.  these  characteristics  result  directly  from 
the  velocities  induced  by  the  strong  tip  vortex.  Although  the  radial 
extant  of  tbe  vortex  sheet  is  depicted  in  Figure  8  as  adding  abruptly,  it 
probably  retain?  soma  compaction  with  tbs  tip  aortas.  Tbs  exact  nature 
of  this  coNMction,  however,  has  bean  difficult  to  distinguish  in  flow 
visualisation  studies .  in  addition  to  Urn  wake  struct***  <  shown  in 
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Figure  8  for  one  blade,  similar  wake  structures  for  other  blades  also 
exist,  with  the  aggregate  forming  the  complete  wake  representation. 

The  schematic  in  Figure  8  is  representative  of  the  three-dimensional 
wake  pattern  which  would  be  observed  if  smoke  were  emitted  from  the  blade 
(i.e.j  in  the  rotating  system).  To  take  advantage  of  the  symmetrical 
nature  of  the  wake  of  a  hovering  rotor  as  well  as  to  facilitate  the 
acquisition  of  quantitative  data,  the  wake  for  this  investigation  was 
observed  by  emitting  smoke  externally  from  the  blades  (i.e.,  in  the 
nonrotating  system).  Smoke  was  emitted  from  smoke  rakes  in  a  single  plane 
and  the  flow  patterns  were  photographed,  as  shown  for  the  two-bladed  rotor 
in  Figure  9«  In  this  manner,  a  two-dimensional  cross  section  of  the  wake 
was  recorded.  The  cross  sections  of  the  tip  vortices  appear  as  circles 
in  which  the  central  regions  are  clear  of  smoke  due  to  the  local  centrif¬ 
ugal  field  which  forces  the  smoke  particles  radially  outward.  The  center 
of  the  circular  cross  sections  are  interpreted  as  the  centers  of  the 
vortex  core.  The  vortex  sheet  cross  sections  sire  indicated  by  the 
discontinuities  present  in  the  smoke  filaments  passing  through  the  inner 
region  of  the  rotor  wake. 

The  photographic  wake  data  were  analyzed  for  selected  test  conditions 
to  determine  the  principal  wake  geometry  characteristics .  The  conditions 
were  selected  so  as  to  permit  assessment  of  the  effects  of  the  primary 
test  parameters.  Radial  and  axial  wake  coordinates  were  determined  from 
the  photographs  as  functions  of  the  wake  azimuth  angle  (’Aw),  which  is 
equivalent  to  the  blade  azimuth  travel  (  <A  =  fit)  from  the  time  it 
generates  the  vortex  cross  section.  For  example,  on  the  right  side  of 
Figure  9>  the  uppermost  tip  vortex  and  vortex  sheet  cross  section  were 
shed  by  blade  2,  which  has  travelled  180  deg  from  the  time  it  passed 
through  the  plane  of  smoke  on  the  right  side  (reference  side).  The 
following  tip  vortex  and  vortex  sheet  cross  section  were  shed  by  blade  1 
the  previous  time  it  passed  through  the  reference  side  of  the  smoke  plane, 
and  thus  the  wake  azimuth  angle  for  these  cross  sections  is  360  deg. 
Likewise,  the  azimuth  angle  of  the  third  cross  section,  shed  by  blade  2, 
is  5^0  deg.  It  should  be  noted  that  the  cross  sections  in  the  visible 
wake  near  the  rotor  remain  approximately  in  the  same  plane  (rotor  wake 
tangential  velocities  are  small).  The  azimuth  angle,  <AW,  was  used  as  the 
third  coordinate  in  the  wake  geometry  analysis.  In  Figure  10,  a  schematic 
of  the  wake  of  Figure  9  and  the  wake  coordinate  system  is  presented.  For 
a  stable,  hovering  rotor  wake,  the  radial  and  axial  coordinates  at  a 
given  azimuth  angle  are  equivalent  for  each  blade  due  to  symmetry.  Thus 
the  complete  coordinate  system  of  the  wake  for  a  given  test  condition  was 
determined  by  the  following  procedure. 
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A  transparent  grid  template  was  constructed  from  the  photograph  of 
the  reference  grid  shown  in  Figure  6.  With  the  grid  template  as  an  over¬ 
lay,  the  radial  and  axial  coordinates  of  the  wake  from  several  blades  in 
a  single  photograph  were  determined  along  with  the  corresponding  wake 
azimuth  coordinates.  This  was  repeated  for  a  sequence  of  photographs 
taken  with  the  rotor  at  a  series  of  prescribed  rotational  positions.  A 
sample  sequence  is  presented  in  Figure  11,  in  which  the  rotor  rotational 
positions  are  designated  by  the  azimuth  position,  ,  of  the  blade  which 
most  recently  passed  through  the  reference  plane  containing  the  smoke . 
Considering  the  known  azimuth  positions  of  the  reference  blade  and 
relating  each  vortex  cross  section  to  the  appropriate  blade,  the  radial 
and  axial  coordinates  were  obtained  from  the  photographs  in  Figure  11  for 
the  following  wake  azimuth  angles  : 

Reference  Blade  Wake  Azimuth, 

Azimuth,  ,  deg  ^w  ,  deg 


0 

0,180,360,540 

15 

15,195,375,555 

30 

30,210,390,570 

6o 

60,240,420,600 

90 

90,270,450,630 

120 

120,300,480,660 

150 

150,330,510 

The  coordinate  results  from  the  series  of  photographs  were  then  plotted 
as  functions  of  ^ w  as  shewn  in  Figures  12  and  13.  To  facilitate  the 
comparison  of  wake  geometries  from  varying  rotors  and  test  conditions,  the 
radial  and  axial  coordinates  were  nondimensionalized  by  the  rotor  radius, 
and  differences  in  axial  coordinates  due  to  blade  coning  were  eliminated 
by  using  the  blade  tip  as  the  reference  (z<j  instead  of  z).  The  range  of 
wake  azimuth  angles  for  which  data  could  be  acquired  was  limited  by  the 
visibility  of  the  smoke.  For  example,  for  two-  and  four-biaded  rotors, 
generally  only  1  to  2  tip  vortex  revolutions  per  blade,  were  visible. 

For  six-  and  eight-bladed  rotors,  less  than  one  revolution  from  each  blade 
was  visible.  However,  it  will  be  shown  that  rotor  performance  is  mainly 
sensitive  to  this  near  wake  geometry  and  insensitive  to  reasonable  extra¬ 
polations  of  this  geometry  made  to  define  the  far  wake.  Thus,  knowledge 
of  the  exact  positioning  of  the  far  wake  elements  is  not  essential  to  the 
objectives  of  the  rotor  performance  portion  of  this  study. 
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Since  the  cross  sections  of  the  vortex  sheet  from  each  blade  are 
essentially  lines  rather  than  discrete  points,  as  is  the  case  for  the  tip 
vortex  cross  sections,  the  following  procedure  for  transforming  the 
photographic  data  for  the  vortex  sheets  to  coordinate  form  was  found  to 
be  convenient.  Assuming  the  vortex  sheet  cross  sections  to  be  linear, 
the  axial  position  of  a  vortex  sheet  at  a  given  azimuth  can  be  defined  by 
two  points.  For  simplicity,  the  two  points  selected  were  the  imaginary 
extensions  of  the  cross  section  to  r  =  0  at  one  end  and  to  r  =  1.0  at  the 
other  end,  as  shown  in  Figures  10  and  13.  These  two  points  establish  the 
intercept  at  the  axis  of  rotation  and  the  slope  of  the  vortex  sheet. 

Using  this  procedure,  it  is  unnecessary  to  reduce,  for  each  smoke  filament, 
the  coordinates  of  their  intercepts  with  the  vortex  sheet  cross  sections 
if  the  following  assumption  is  made.  It  was  assumed  that  the  radial 
position  of  such  intercepts  (e.g.,  point  A  of  Figure  10)  is  linearly 
proportional  to  the  radial  coordinate  of  the  intersection  of  the  vortex 
sheet  with  the  vortex  sheet  boundary  (point  B).  With  the  exception  of 
the  immediate  vicinity  of  the  blade  (<^w  <  2 7r/b ) ,  the  vortex  sheet  boundary 

was  assumed  to  be  equivalent  to  the  boundary  formed  by  the  locus  centers 

of  the  tip  vortex  cross  sections  (equivalently,  tip  vortex  streamline). 

For  'fry  less  than  the  blade  spacing,  27r/b,  the  boundary  was  faired  from 

the  point  of  maximum  circulation  on  the  blade  (point  D  of  Figure  10)  to 

the  tip  vortex  boundary.  The  constant  of  proportionality  was  assumed  to 
be  the  ratio  of  the  radial  position  of  the  origination  of  the  smoke 
filament  streamline  at  the  blade  (point  C)  to  the  radial  position  of  the 
vortex  sheet  boundary  at  the  blade  (point  T5).  That  is, 


It  was  found  that  this  was  a  reasonable  approximation  for  use  in  the 
theoretical  wake  model  for  rotor  performance  calculations  to  be  presented 
in  a  later  section. 


TEST  DATA  ACCURACY 


Static  data  repeatability  for  thrust  and  torque  was  determined  from 
repeated  calibrations  of  the  strain  gages  made  while  determining  the 
calibration  derivatives  discussed  in  the  Test  Procedures  section.  The 
repeatability  values,  represented  by  two  standard  derivations,  are 
listed  below : 
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Thrust,  lb  Torque,  ft- lb 

Static  Repeatability:  ±0.0087  -0,1130 

The  dynamic  data  repeatability  was  established  by  considering  the  range 
of  Ct/ct  and  C q/cr  measurements  observed  from  consecutive  test  points  at 
a  given  tip  speed  for  each  test  configuration.  Normally,  four  test 
points  were  available  from  two  consecutive  runs  between  which  the  rotor  was 
stopped  and  restarted.  The  mean  range  was  established  and  related  to  the 
standard  deviation  using  the  procedures  outlined  by  Hoel  (Reference  4, 
p.  241).  The  results  obtained  are  given  below. 


b=2, 

b=8 

Average 

fiR=525  fps 

HR =700  fps 

Dynamic  Data  Repeatability,  Cj/<r 

±0. 00024 

-0.00044 

±0. 00015 

cQ/o- 

±0.000022 

±0.000045 

± 0.000012 

Ry  choosing  consecutive  test  points  in  the  above  analysis,  variations  due 
to  differences  in  collective  pitch  setting  were  eliminated.  These 
differences  were  significant,  as  will  be  shown  in  the  presentation  of  the 
performance  data.  However,  they  did  not  significantly  influence  the  rotor 
thrust  -  torque  relationship. 

The  estimated  accuracies  with  which  the  parameters  determining  a 
given  test  condition  could  be  set  are  as  follows : 


Parameter 


Accuracy 


Collective  Pitch,  6j 5 
Tip  Speed,  Hr 
Rotor  Height,  Zq/r 


-0.2  deg 
*1  fps 
±0.03  R 


The  estimated  accuracy  and  repeatability  of  the  tip  vortex 
coordinates  are  listed  below. 


Wake  Coordinate 

Accuracy 

Repeatability 

Azimuth, {IJV 

-3  deg 

— 

Radial,  f 

±0.005 

-  ±0.01 

Axial,  zT 

±0.005 

~  ±0.01 

13 


The  accuracy  of  the  radial  and  axial  coordinate  represents  the  degree 
of  accuracy  to  which  coordinates  of  a  wake  point  could  be  measured  from 
the  grid  system  used.  The  repeatability  represents  the  normal  repeat¬ 
ability  of  a  wake  point  in  a  series  of  photographs  taken  during  a  single 
test  condition  (e.g.,  see  Figure  12).  This  repeatability  pertains  only 
to  the  tip  vortex  in  the  near  wake  region  of  the  rotor  (to  be  discussed 
later'.  Figure  14  shows  the  repeatability  of  the  tip  vortex  coordinates 
for  a  test  condition  as  repeated  on  three  different  dates. 


DISCUSSION  OF  EXPERIMENTAL  ROTOR  PERFORMANCE  RESULTS 


The  rotor  performance  test  data  were  transformed  to  graphical  form 
for  the  various  combinations  of  number  of  blades  (solidity),  blade  twist, 
aspect  ratio,  tip  speed,  and  simulated  ground  height.  The  resulting 
graphs  were  analyzed  to  assess  the  separate  influence  of  each  of  these 
parameters  on  model  rotor  hover  performance.  The  performance  character¬ 
istics  of  the  model  rotors  are  presented  in  nondimens ional  form  in 
Figures  15  through  18  for  each  combination  of  blade  twist  and  aspect 
ratio  tested.  This  series  of  graphs,  which  contains  data  for  out-of¬ 
ground-effect  operation  (Zq/r  =  3.5),  also  includes  results  for  various 
tip  speeds. 

Thrust  (Collective  Pitch)  Limits 

As  noted  previously,  maximum  thrust  (collective  pitch)  was  limited, 
for  each  rotor  configuration,  by  the  rapid  increase  in  noise  level  when 
a  specific  tip  speed  was  exceeded.  Oscillograph  records,  displaying  a 
signal  produced  by  torsional  strain  gages  mounted  at  the  root  end  of  a 
blade,  were  analyzed  to  determine  the  frequency  of  the  blade  torsional 
response.  The  frequency  ( oj  =  approximately  12  cycles  per  rotor  revolu¬ 
tion)  was  found  to  agree  with  both  the  measured  acoustic  frequency  and  an 
estimate  of  the  first  natural  frequency  of  the  blade  in  torsion  when 
nondimensionalized  by  fl.  Thus,  it  was  concluded  that  the  performance 
boundary  had  the  characteristics  of  incipient  stall  flutter.  Previous 
stall  flutter  investigations  (e.g.,  Reference  5)  have  shown  that  the 
thrust  boundary  (hereafter  referred  to  as  the  stall  flutter  boundary) 
is  lowered  with  an  increase  in  the  following  parameter  which  is  the 
inverse  of  the  reduced  frequency  parameter : 

STALL  FLUTTER  PARAMETER  =  4~- 

c  c o 
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The  test  results  indicated  that  the  stall  flutter  boundaries  were  lowered 
with  increased  tip  speed, fiR,  and  decreased  chord,  c,  indicating  apparent 
agreement  with  the  above  relation.  The  stall  flutter  boundaries  were 
also  lowered  with  decreasing  blade  twist  due  to  the  earlier  stall  of  the 
blade  tip  section. 

Effect  of  Collective  Pitch 


The  results  of  Figures  15  to  18  are  presented  in  terms  of  rotor 
thrust  versus  rotor  torque  to  minimize  scatter  introduced  by  the  accuracy 
with  which  collective  pitch,  6j 5,  could  be  set.  Nominal  values  of 
constant  collective  pitch  are,  however,  indicated  by  the  dashed  lines  in 
each  figure.  In  addition,  representative  variations  of  C^/cr  and  Cq/o-  with 
9y<j  are  presented  in  Figure  19  corresponding  to  the  faired  data  presented 
in  Figure  17.  As  anticipated,  rotor  torque  increases  rapidly  at  the 
higher  collective  pitch  levels  due  to  divergence  of  the  airfoil  drag 
characteristics.  However,  no  corresponding  fall-off  in  thrust,  an 
indication  of  lift  stall,  is  evident  for  the  range  of  pitch  values  tested. 

Effect  of  Solidity  and  Number  of  Blades 

Rotor  solidity  was  changed,  for  each  blade  design,  by  varying  the 
number  of  blades  from  2  to  8  as  indicated  in  Figures  15  through  18.  Since 
each  of  these  figures  represents  rotors  with  fixed  blade  chord,  radius, 
and  tip  speed,  Ct/ct  and  Cq/ct  are  directly  representative  of  thrust  and 
torque  per  blade  (blade  thrust  and  torque  loading).  Comparison  of  the 
results  for  varying  numbers  of  blades  within  each  figure  demonstrates  the 
improved  blade  efficiency  (thrust/blade  (C?  /<r)  per  torque  /blade  (Cq/o-)) 
with  decreasing  solidity.  The  improved  blade  efficiency  occurs  for  two 
reasons.  First,  and  most  important,  at  a  given  blade  loading,  the  fewer 
the  number  of  blades  (or  the  smaller  the  solidity),  the  lower  will  be  the 
total  thrust  and  disc  loading  of  the  rotor.  As  a  result,  the  average 
downwash  (which  is  a  measure  of  energy  expended)  induced  by  the  rotor 
wake  will  be  lower  and  hence  the  blade  induced  drag  will  be  less.  A 
second  smaller  effect  arises  from  the  reduced  local  interference  caused 
by  the  tip  vortex  shed  by  one  blade  on  the  loading  of  the  following  blade. 
This  will  be  discussed  further  in  a  later  section  of  this  report. 

Although  lower  solidity  at  constant  aspect  ratio  implies  greater 
blade  efficiency,  there  is  a  limit  to  which  it  also  implies  greater  rotor 
efficiency.  This  is  shown  in  Figure  20,  in  which  the  total  rotor  perfor¬ 
mance  coefficients  (cT  and  Cq)  are  presented  for  the  blade  design  and  tip 
speed  corresponding  to  Figure  16(a).  Since  the  blade  radius  and  tip  speed 


15 


are  fixed  in  Figure  20,  and  Cq  are  representative  of  total  rotor  thrust 
and  torque.  Figure  20  shows  that  the  rotor  efficiency  (rotor  thrust  per 
rotor  torque)  is  improved  with  lower  solidity  only  up  to  a  point  which  is 
determined  by  blade  stall  considerations. 

Effect  of  Aspect  Ratio  at  Constant  Solidity 

Of  particular  interest  to  the  rotor  designer  is  the  trade-off  in 
performance  between  chord  and  number  of  blades  while  maintaining  ccr.3tant 
rotor  solidity  (total  blade  area  and  disc  area  held  constant).  The 
experimental  results  comparing  the  hover  performance  for  eight  18.2- 
aspect-ratio  blades  (c  =  1.47  in.)  and  six  13.6-aspect-ratio  blades 
(c  =  1.96  in.)  at  a  constant  solidity  of  0.140  are  presented  in  Figure  21. 
Over  the  thrust  range  tested  (i.e.,  up  to  the  stall  flutter  boundary), 
the  results  are  essentially  equivalent  for  the  two  configurations.  The 
existence  of  the  stall  flutter  boundary  prohibited  the  investigation  of 
the  trade-off  of  number  of  blades  and  chord  at  conditions  associated  with 
deep  penetration  into  stall.  The  eight  narrow  chord-blades  exhibited 
stall  flutter  at  lower  performance  levels  than  the  six  wide-chord  blades. 
This  implies  that  the  aeroelastic,  rather  than  the  aerodynamic,  character¬ 
istics  of  the  blades  may  ultimately  be  the  determining  factor  in  selecting 
blade  aspect  ratio. 

Effect  of  Blade  Twist 


The  effect  of  blade  linear  twist  rate  on  experimental  model  rotor 
performance  is  shown  in  Figure  22  for  two-  and  six-bladed  rotors  operating 
at  a  tip  speed  of  700  fps.  For  clarity,  only  the  faired  curves  through 
the  experimental  data  from  Figures  15  to  17  are  presented.  An  improvement 
in  performance  is  noted  at  the  higher  thrust  levels  as  blade  twist  is 
increased.  For  the  six-bladed  rotor  most  of  the  improvement  is  due  to 
the  initial  8  degrees  of  twist.  For  the  two-bladed,  -l6-degree-twist 
rotor  blades,  the  higher  profile  drag  at  low  thrust  and  the  improved 
performance  in  the  stall  region  result  in  a  crossover  in  blade  efficiency 
relative  to  the  blades  of  lower  twist. 

Effect  of  Tip  Speed 

The  effect  of  variations  in  tip  speed  (from  525  to  700  fps)  on  rotor 
performance  is  shown  in  Figures  23  and  24.  As  in  Figure  22,  only  the 
faired  experimental  curves  are  presented.  Performance  results  are  pre¬ 
sented  for  two-  and  six-bladed  rotors  with  0  and  -l6-degree-twist  blades. 
The  influence  of  compressibility,  which  results  in  decreasing  performance 


16 


with  increasing  tip  speed,  is  evident  for  all  rotors.  The  influence  is 
less  for  the  -l6-degree-twist  blades  because  the  high-Mach-number  tip 
sections  are  operating  at  lower  angles  of  attack  relative  to  the  tip 
sections  of  the  zero-twist  blades. 

Ground  Effect 

The  relative  position  of  the  rotor  with  respect  to  the  ground  was 
varied  by  raising  or  lowering  the  ground  plane.  The  rotor  configuration 
and  collective  pitch  setting  were  held  fixed  during  a  series  of  simulated 
rotor  height  variations.  The  effect  of  rotor  height  above  ground  (ground 
effect)  on  rotor  performance  coefficients  is  shown  in  Figure  25  for  six- 
and  eight-bladed  rotors  with  a  -8-degree  twist  and  an  aspect  ratio  of 
18.2.  Data  are  presented  for  variations  in  rotor  height  from  3.5  R,  which 
is  essentially  out  of  ground  effect  (OGE),  to  0.6?  R,  which  is  well  in 
ground  effect  (IGE).  The  OGE  data  differs  slightly  from  that  presented 
previously  because  the  plastic  tip  blades  were  used  for  this  series  of 
data;  however,  this  should  have  no  effect  on  the  relative  performance 
between  the  1(3!  and  OGE  conditions.  It  is  shown  in  Figure  25  that  the 
variation  of  thrust  with  torque  is  essentially  linear  as  rotor  height  is 
varied.  Thrust  augmentation,  which  is  defined  as  the  increase  in  rotor 
thrust  at  constant  torque  for  1(21  operation  over  that  for  OGE  operation, 
was  obtained  using  the  following  relation: 


(Ct/o)I6 

(Ct/<t)c 


°®EJ  CONS  i  ANT  Cq/0 r 


The  thrust  augmentation  results  are  presented  in  Figure  26,  where  it  is 
shown  that  the  effect  of  rotor  height  is  virtually  independent  of  both 
variations  in  solidity  corresponding  to  six  and  eight  blades  and 
variations  in  OGE  thrust  level.  The  thrust  at  the  minimum  rotor  height 
tested  was  increased  by  approximately  18^  over  the  thrust  for  OGE  opera¬ 
tion.  It  should  be  noted  that  the  results  presented  are  for  -8-degree 
linear  twist  blades.  The  detailed  investigation  of  ground  effect  was  not 
a  primary  objective  of  this  study  and  thus  was  not  investigated  for  all 
rotor  configurations.  Results  of  recent  tests  at  UARL  and  Sikorsky 
Aircraft  indicate  a  variation  in  thrust  augmentation  with  blade  twist. 

In  addition  to  the  variation  in  rotor  height  discussed  above,  the  -8- 
degree-tvist  blades  (AR  »  18.2)  were  tested  I®  in  a  manner  to  simulate 
typical  full-scale  whirl-stand  operation.  For  these  conditions,  the  rotor 


was  1.67  R  above  the  ground  plane  and  was  tested  with  and  without  the 
presence  of  a  model  whirl  3tand  (IGE  and  IGE/WS).  The  following  tabula¬ 
tion  consists  of  representative  averages  of  the  thrust  augmentation 
results. 

Thrust  Augmentation  Ratio,  Tjqjj/Toqe 


Tip  Speed 

No.  of 

(fps) 

Blades 

IGE 

IGE/WS 

700 

2,  6,  8 

1.03 

1.035 

525 

2,  6,  8 

1.02 

1.025 

Although  the  scatter  in  the  individual  Tjqjj/Tqqe  results  used  in  determin¬ 
ing  the  above  averages  was  approximately  ±l£,  taken  in  aggregate,  the 
data  indicated  a  slight  increase  in  thrust  augmentation  due  to  the 
presence  of  the  whirl  stand. 


DISCUSSION  OF  EXPERIMENTAL  ROTOR  WAKE  GECMETRY  RESULTS 


Sample  Photographs 

In  addition  to  the  sequence  of  photographs  previously  presented  in 
Figure  11  for  a  two-bladed  rotor,  similar  sequences  for  four*  six,  and 
eight  blades  are  presented  in  Figures  27  through  29.  These  photographs, 
for  the  zero  twist  rotor  operating  at  a  tip  speed  of  700  fits  and  a 
collective  pitch  setting  of  8  degrees,  show  the  time  histcry  of  the  wake 
as  the  reference  blade  rotates  to  various  azimuth  positions.  The  tip 
vortex  cross  sections  in  the  "near  wake"  (portion  of  wake  within  one  or 
two  wake  revolutions  beneath  the  rotor)  are  clearly  evident,  as  are  the 
discontinuities  in  the  inboard  smoke  filaments  which  identify  the  inboard 
vortex  sheet  locations.  Also  evident  are  the  rapid  contraction  of  the 
wake  and  the  decreased  axial  wake  spacing  with  increasing  number  of 
blades. 

A  characteristic  of  ail  the  photographs  taken  during  this  investiga¬ 
tion  is  the  absence  of  the  tip  vortex  cross  sections  beyond  those 
corresponding  to  approximately  three  to  four  blade  passages.  That  is,  no 
more  than  four  vortex  cross  sections  normally  appear  clearly  cm  any 
given  photograph.  Thus,  the  limiting  azimuth  for  visualization  of  the  tip 
vortex  varied  with  number  of  blades  according  to  the  following: 
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No.  of  Blades 

Azimuth  Spacing 
Between  Blades 
(deg) 

Approximate  Azimuth  Limit 
(Azimuth  Spacing  x  4) 
(deg) 

2 

180 

720 

4 

90 

360 

6 

60 

240 

8 

45 

180 

At  first  it  was  believed  that  the  tip  vortex  cross  sections  became 
indistinguishable  due  to  diffusion  of  the  smoke  by  turbulent  wiring  as 
mentioned  by  Gilmore  in  Reference  6.  However,  close  scrutiny  of  the 
photographs  and  high-frame-speed  movies  resulted  in  the  conclusion  that 
the  far-wake  region  of  a  hovering  rotor  is  unstable  or,  at  best,  neutrally 
stable.  Whether  the  tip  vortex  undergoes  a  form  of  viscous  dissipation 
(decay)  or  vortex  breakdown  (bursting)  as  characterized  by  certain  fixed 
wing  tip  vortices  (Reference  7)  is  conjecture  at  this  time.  However,  the 
results  of  this  investigation  indicate  a  definite  departure  from  the 
classical  concept  of  a  smoothly  contracting  wake  below  the  rotor. 

Evidence  of  this  is  shown  in  Figure  11,  in  which  the  fourth  vortex  cross 
section  proceeds  to  travel  radially  outward  (note  particularly  the 
photographs  for  ^=90  and  120  deg)  until  it  is  no  longer  visible  at 
♦  a  1^0  deg.  Perhaps  a  dramatic  example  of  the  instability  is 
presented  in  the  series  of  Figo.3  30,  in  which  the  fourth  vortex  cross 
section  moves  radially  out  and  eventually  to  the  aide  of  the  third. 
Although  it  la  recognized  that  a  small  perturbation  such  as  a  small 
amount  of  ambient  wind  or  a  slight  blade-out-of- track  may  be  necessary  to 
precipitate  the  instability,  this  la  believed  to  be  an  academic  consider* 
atioo  since  full-scale  rotcar  operation  is  certainly  subject  to  greater 
disturbances  than  those  present  under  the  laboratory  conditions  of  this 
Investigation.  More  photographic  examples  and  a  more  complete  discussion 
of  the  unstable  nature  of  the  far  wake  will  be  presented  in  a  later 
section  of  this  report. 

Sample  enlarged  photographs  for  combinations  of  rotor  configurations 
and  test  conditions  representative  of  vsriations  of  the 
primary  teat  parameters  are  presented  in  Figures  31  through  4o.  Three 
photographs  for  IGE  and  whirl-stand  conditions  are  also  included.  Single 
photographs  for  each  configuration  and  condition  are  shown  with  the 
reference  blade  at  the  zero  azimuth  position.  Note  that  the  top  Up 
vortex  cross  section  in  each  photograph  corresponds  to  the  tip  vortex  shed 
from  the  previous  blade.  A  time  exposure  (Figure  32)  of  the  wake  is 
Included  to  shew  the  streamline  patterns  of  the  smoke  filaments.  Note  that 
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this  photograph  is  for  the  same  condition  as  one  of  the  previous  stop- 
action  photographs  presented  in  Figure  31*  Sample  photographs  for  the 
rotor  in  a  simulated  whirl-stand  environment  and  IGB  are  presented  in 
Figures  38  through  40.  The  radial  expansion  of  the  wake  as  it  approaches 
the  ground  plane  is  evident  in  these  latter  photographs.  One  of  the 
observations  made  during  review  of  these  and  other  photographic  data  was 
that  the  tip  vortex  never  passed  above  the  blades.*  Ffcotographic 
evidence  of  the  location  of  the  tip  vortex  of  a  full-scale  rotor  above  the 
plane  of  rotation  was  shown  in  Reference  1,  where  it  was  attributed  to  a 
3 mall  amount  of  wind.  The  absence  of  this  phenomenon  in  the  model  tests 
is  ascribed  to  the  controlled  laboratory  ccndi tions.  The  wake  character¬ 
istics  observed  in  the  flow  visualization  photographs  are  discussed  in 
more  detail  in  the  following  section. 

Tip  Vortex  Coordinates 


Sample  results  obtained  from  the  reduction  of  the  flow  visualisation 
photographs  to  coordinate  farm,  are  presented  in  Figures  41  through  52 
far  the  tip  vertex.  Note  that  in  order  to  present  the  axial  and  radial 
coordinates  on  the  same  scale,  and  since  the  experimental  tip  vortices 
were  always  located  below  the  rotor  tip  path  plana,  the  negative  of  the 
axial  coordinate  vitn  respect  to  the  bl  tde  tip  (-1^)  was  plotted.  The 
symbols  at  e  specific  azimuth  represent  an  average  of  ell  the  data 
reduced  from  multiple  photographs  for  a  given  rotor  operating  condition. 
The  repeatability  of  the  redial  and  axial  coordinates  was  generally 
within  ^O.Cl  R  la  the  stable  neer-wakt  region.  The  flight  conditions  for 
which  data  are  presented  were  selected  from  over  70  analysed  conditions 
to  provide  typical  results  indicating  the  effect  of  independent  variations 
of  the  primary  teat  peremeters.  Since  the  data  in  Figures  4i  through  $2 
correspond  to  specific  test  points,  it  is  suggested  that  they  be  used 
only  to  exea$lify  trends  but  not  for  final  wake  generalisation  purposes. 
Generalized  plots  based  on  data  from  all  reduced  te*t  conditions  will  be 
presented  in  a  later  section. 


•In  some  of  the  photographs  showing  the  time  history  of  the  wake  {e.g., 
Figures  27  and  28),  e  tip  «arte*  cross  section  appears  to  be  above  the 
blade.  Ibis  is  attributed  to  the  viewing  angle  of  the  camera  which  waa 
placed  beneath  the  rotor  plane.  Since  the  blades  in  these  photographs 
move  out  of  the  plane  of  the  smoke,  their  position  is  distorted  relative 
to  the  wake  cross  sections  in  that  plane. 
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General  Wake  Features 


Several  general  features  of  the  tip  vortex  geometry  are  evident  in 
Figures  4l  through  52.  When  an  element  of  the  tip  vortex  is  shed  from  a 
blade,  its  rate  of  axial  displacement  is  very  low  until  it  passes  beneath 
the  following  blade  (at  -  360  deg/b).  At  that  point,  the  tip  vortex 
clement  lies  radially  inboard  of  the  tip  vortex  of  the  following  blade 
and  thus  experiences  a  large  downward  induced  velocity  from  that  vortex. 
The  axial  transport  velocities  before  and  after  the  passage  of  the 
following  blade  are  fairly  constant  in  the  near  wake,  as  can  be  seen  from 
the  substantially  linear  variations  of  the  axial  displacement,  z-p,  with 
wake  azimuth  angle  in  these  regions.  The  radial  displacement,  r,  of  the 
tip  vortex  decays  in  an  apparently  exponential  manner  as  the  wake  azimuth 
is  increased.  The  fairings  of  the  data  in  the  figures  are  based  on  these 
general  wake  characteristics . 

Effect  of  Thrust  Level  (Collective  Pitch)  for  Fixed  Number  of  Blades 

The  tip  vortex  coordinates  for  a  two-bladed  rotor  operating  at  three 
thrust  levels  corresponding  to  three  collective  pitch  values  are  shown  in 
Figure  4l.  Figure  42  shows  the  coordinates  for  a  six-bladed  rotor  for  a 
similar  variation.  The  axial  displacement  and  radial  contraction  of  the 
tip  vortex  is  shown  to  increase  with  increasing  thrust  level.  It  should 
be  noted  that  this  result  differs  from  the  generalized  wake  model 
assumed  by  Rorke  and  Wells  (Reference  8),  which  was  based  on  limited 
experimental  data. 

Effect  of  Number  of  Blades  at  Constant  Blade  Loading 

The  effect  on  tip  vortex  geometry  of  varying  the  number  of  blades 
while  holding  the  blade  loading  (or  Cip/oO  constant  is  shown  in  Figure  43 
for  the  zero-degree-twist  rotor.  The  rates  of  axial  displacement  for 
the  two-  and  six-bladed  rotors  are  equal  up  to  the  azimuth  position 
corresponding  to  that  of  the  following  blade.  Beyond  this  point  the  tip 
vortex  of  the  six-bladed  rotor,  which  has  a  greater  total  thrust,  travels 
downward  at  a  faster  rate.  The  radial  contraction  for  the  six-bladed 
rotor  is  also  greater  than  that  of  the  two-bladed  rotor.  Similar  results 
are  shown  in  Figure  44  for  a  -8-degree-twist  rotor. 
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Effect  of  Number  of  Blades  at  Constant  Dice  Loading 


The  effect  on  tip  vortex  geometry  of  varying  the  number  of  blades 
while  holding  the  rotor  disc  loading  (or  CT)  constant  is  shown  in  Figure 
45.  The  rate  of  axial  displacement  is  greater  for  the  two-bladed  rotor, 
which  has  the  greater  blade  loading,  up  to  the  azimuth  position  corre¬ 
sponding  to  the  following  blade.  Beyond  this  point  the  displacement  rates 
are  equivalent.  Also,  the  radial  coordinates  are  equivalent  at  constant 
disc  loading. 

Effect  of  Blade  Twist 

To  show  the  independent  effect  of  blade  twist  rate  on  tip^firtex 
geometry,  the  tip  vortex  coordinates  ha^e  been  plotted  in  Figure?j46  for 
a  two-bladed  rotor  at  a  constant  thrust  level.  Although  the  radial 
coordinates  are  independent  of  twist  variations,  the  axial  dispBcement 
decreases  with  increasing  twist  rate  --  particularly  between  -8^nd  -16 
degrees.  Similar  results  are  shown  in  Figure  47  for  a  six-bladed  rotor. 

Effect  of  Aspect  Ratio  at  Constant  Solidity 

The  effect  of  aspect  ratio  on  tip  vortex  geometry  is  shown  in 
Figure  48,  in  which  the  data  “or  six  wide-chord  blades  are  compared  at  a 
constant  thrust  level  with  the  data  for  eight  narrow-chord  blades  having 
the  same  solidity  ratio.  The  radial  coordinates  and  the  axial  displace¬ 
ment  rates  are  approximately  equivalent.  The  axial  displacement  after 
the  following  blade  passes  is  greater  for  the  high-aspect-rat,io  rotor 
due  +0  che  lesser  blade  spacing  of  the  eight-bladed  rotor. 

Effec  C  v>A  •  Tip  Speed 

Tip  vortex  geometry  is  shown  in  Figures  49  and  50  to  be  independent 
of  tip  speed.  Data  for  two-  and  six-bladed  rotors  are  presented  for  tip 
speeds  of  525 »  600,  and  ^CO  fps,  which,  as  previously  mentioned, 
correspond  to  the  tip  Mach  number  range  of  0.46  to  0.6l. 

Effect  of  Ground  and  Whirl  Stand 


The  effect  of  the  simulated  ground  and  whirl  stand  on  tip  vortex 
geometry  is  shown  in  Figures  51  and  52  for  two-  and  six-bladed  rotors. 
For  each  rotor,  the  collective  pitch  setting  was  held  constant,  and  the 
ground  plane  position  was  changed  from  the  OGE  to  the  IGE  position  and 
the  whirl-stand  model  inserted.  As  shown  in  these  figures,  the  axial 
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displacement  of  the  tip  vortex  is  not  sensitive  to  the  whirl  stand  and 
corresponding  ground  plane  position  (Zq/r  =  1.67).  Conflicting  results 
were  obtained  for  the  radial  coordinates  of  the  two-  and  six-blade  tip 
vortices.  For  the  two-bladed  rotor,  the  presence  of  the  ground  or  ground 
and  whirl  stand  resulted  in  an  expansion  of  the  tip  vortex  boundary 
(Figure  51).  However,  no  difference  was  evident  for  the  six-bladed  rotor 
(Figure  52).  It  is  noted  that  these  IGE  results  are  based  on  a  limited 
amount  of  data,  and  further  investigation  is  recommended.  Also,  the  rotor 
height  above  ground  for  the  simulated  whirl-stand  condition  does  not 
correspond  to  what  would  be  considered  a  severe  ground  effect  condition 
(see- Figure  26).  Investigation  of  the  effect  of  ground  height  on  wake 
geometry  was  beyond  the  scope  of  this  program. 

Gerer^fc^pd  Tip  Vortex  Geometry  Results 

It  hks  been  shown  in  the  preceding  section  that  the  tip  vortex  flow 
visualization  data  can  be  well  represented  by  a  series  of  straight  lines 
for  the  axial  coordinate  and  by  an  exponential  function  for  the  radial 
coordinate.  The  parameters  defining  these  curves  were  determined  for  all 
conditions  that  were  reduced  and  cross  plotted  in  an  attempt  to  generalize 
the  tip  vortex  geometry  in  terms  of  fundamental  rotor  parameters.  The 
procedures  followed  and  results  obtained  are  described  below.  It  is 
recognized  that  some  of  the  conclusions  reached  in  interpreting  these  data 
depend  on  the  manner  in  which  the  data  are  faired.  The  fairings  used 
represent  our  interpretation  of  the  data  and  are  believed  to  be  consistent 
and  reasonable,  taking  into  consideration,  as  they  do,  all  of  the  data 
.ilabe  for  the  different  blade  designs. 


Vortex  Axial  Coordinate 


All  axial  coordinate  data,  including  those  of  Figures  4l  through  52, 
have  been  approximated  by  the  following  relations: 

k|  'I'w  for  0  <  '/'w  <  ^ 

(2) 

+  k2^w  ~  for 
The  constants  k^  and  k2  are  defined,  as  follows : 

k|  orl<2=  -mr'°n'*’w 

2tt 

for  k,  :  A<frw  =  ~y 
for  k2  :  A<//w  "  «//w 
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kl  applies  prior  to  the  azimuth  position  of  the  following  blade  (27r/b)  and 
kg  applies  thereafter.  It  is  shown  in  the  above  equation  that  k^  and  kg 
also  represent  (l)  the  axial  transport  velocity  of  the  tip  vortex,  vz, 
nondimens ionalized  by  the  rotor  tip  speed  and  (2)  the  tangent  of  the  local 
pitch  angle  of  the  contracted  helix  of  the  tip  vortex,  k^  and  kg 
values  were  determined  by  linear  fairings  through  the  axial  coordinate 
data  as  previously  shown  in  Figures  4l  through  52. 

It  was  indicated  in  Figures  42  and  43  that  for  a  given  blade  design, 
k]_  may  be  uniquely  determined  by  blade  loading  (or  more  precisely,  C-p /or). 
Thus,  the  variation  cf  ki  with  the  parameter  Cf/cr  was  plotted  and  is 
presented  in  Figure  53.  Data  for  all  numbers  of  blades  and  tip  speeds  of 
the  test  are  included  for  each  of  the  four  blade  designs  (twist,  aspect 
ratio  combinations).  The  fairings  of  the  data  which  eliminate  differences 
due  to  number  of  blades  and  tip  speed  are  believed  to  be  justified  in 
light  of  the  experimental  accuracy  of  the  flow  visualization  data.  If 
the  repeatability  of  AzT  in  the  kp  region  and  for  a  given  test  condition 
is  -0.015  (see  Figure  14),  the  following  accuracy  of  kp  results  from 
Equation  (2)  for  rotors  of  varying  numbers  of  blades: 

No.  of  Blades  Accuracy  of  kp 


2  ±0.005 

4  ±0.010 

6  ±0.014 

8  ±0.019 

Essentially  all  of  the  data  points  in  Figure  53  lie  within  the  accuracy 
range  of  the  faired  lines.  It  is  possible  that  the  consistency  of  the 
separation  of  the  data  for  six  and  eight  blades,  -8-degree  twist,  and 
an  aspect  ratio  of  18.2  from  the  faired  line  may  be  attributable  to  the 
use  of  the  plastic  tip  blades  for  these  rotor  configurations.  Based  on 
the  faired  lines,  the  following  is  concluded  regarding  the  kp  wake 
parameters  for  the  tip  vortex: 

1.  k^  increases  linearly  with  blade  loading  (or  Cj/cr), 

2.  At  a  given  CT/cr,  k^  is  independent  of  tip  speed. 

3.  At  a  given  CT/cr,  k^  is  independent  of  number  of  blades  -  at 

least  within  the  range  of  experimental  accuracy  of  this 
investigation. 


?4 


4. 


The  rate  of  k^  with  respect  to  C «j»/<r  is  independent  of  blade 
twist  and  aspect  ratio  (slopes  of  k^  vs  Ct /<t  are  essentially 
constant  in  Figure  53)* 

5.  At  a  given  Ct/ct,  k^  increases  with  blade  twist  rate. 

6.  At  a  given  C t/o-,  k^  is  independent  of  aspect  ratio. 

It  is  noted  that  conclusions  5  and  6  are  qualified  by  the  fact  that  the 
variations  of  the  fairings  are  within  the  experimental  accuracy. 

It  was  indicated  in  Figure  44  that  for  a  given  blade  design,  may 
be  uniquely  determined  by  rotor  thrust  coefficient.  Thus,  the  variation 
of  k2  with  the  parameter  Ct  was  determined,  and  the  results  are  presented 
in  Figures  54  through  57-  Data  for  all  numbers  of  blades  and  tip  speeds 
of  the  test  are  included.  Considering  the  experimental  accuracy,  it  was 
again  possible,  for  each  blade  design,  to  determine  a  single  fairing 
which  eliminates  differences  due  to  number  of  blades  and  tip  speed.  It 
was  found  that  the  k2  data  could  be  faired  in  direct  proportion  to  the 
yc-r/2.  This  enables  the  axial  displacement  velocity  of  the  tip  vortex 
in  the  k2  region  to  be  directly  related  to  the  momentum  inflow  velocity 
as  follows : 


vz  =  D.R  k2  =  ftR  C  y/Cj/Z  =  CVj 


MOM. 


The  proportionality  constant,  C,  varied  from  -1.4l  to  -1.19  with  increas¬ 
ing  twist  (0  to  -16  deg), and  was  not  affected  by  variations  in  aspect 
ratio.  Based  on  the  faired  curves,  the  following  is  concluded  regarding 
the  k2  wake  parameter  for  the  tip  vortex  in  the  stable  wake  region: 

1.  k2  is  linearly  proportional  to^jc^/2.  This  implies  that  the 
axial  displacement  velocity  of  the  tip  vortex  is  linearly 
proportional  to  the  momentum  inflow  velocity. 

2.  For  a  fixed  CT  (or  disc  loading),  k2  is  independent  of  tip 
speed  and  number  of  blades. 

3.  The  constant  of  porportionality  between  k2  and  the “\j CT/2  (or 
between  vz  and  vZMa.)  decreases  with  increasing  blade  twist 
but  is  independents  blade  aspect  ratio,  being  given  by  the 
expression  C  *  -1.4l  -  0.0141  0^  (deg). 
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Tip  Vortex  Radial  Coordinate 


It  was  indicated  in  Figure  44  that  the  tip  vortex  radial  coordinate, 
r,  may  be  uniquely  determined  by  rotor  thrust  coefficient.  Thus  the 
variation  of  the  radial  coordinates  with  CT  at  selected  azimuth  positions 
was  determined,  and  the  results  are  presented  in  Figures  58  through  6l. 

The  data  symbols  in  these  figures  represent  the  mean  of  all  available  data 
points  reduced  at  a  given  azimuth  for  a  specific  rotor  configuration  and 
test  condition.  Data  for  all  numbers  of  blades  and  tip  speeds  are 
included,  and,  again,  it  was  possible  to  represent  the  data  with  faired 
lines  which  eliminate  differences  due  to  number  of  blades  and  tip  speed. 
The  faired  lines  shown  for  a  given  azimuth  are  identical  for  each  blade 
design,  indicating  that  the  radial  coordinates  are  also  insensitive  to 
twist  and  aspect  ratio  variations  at  a  constant  Ct.  For  a  given  azimuth, 
the  radial  coordinate  decreases  linearly  with  increasing  Ct.  Further 
analysis  of  the  faired  lines  shown  in  Figures  58  through  6l  indicated  that 
they  could  be  accurately  represented  by  the  following  equation: 


r  =  A  +  (l-A)e'X^w 


(3) 


provided  that  the  constant  A  was  selected  as  O.78  and  that  the  constant  A 
was  determined  from  the  function  of  rotor  thrust  coefficient  given  in 
Figure  62.  The  values  of  A  and  A  were  selected  to  fit  the  near  wake 
radial  coordinates. 

In  summary,  then,  the  following  is  concluded  regarding  the  radial 
coordinates  of  the  tip  vortex  in  the  stable  wake  region: 

1.  For  a  given  azimuth,  the  radial  coordinates  are  linearly 
proportional  to  disc  loading  (or  CT)  and  decrease  with 
increasing  disc  loading. 

2.  For  a  fixed  disc  loading  (or  Cx),  the  radial  coordinates  are 
independent  of  tip  speed,  number  of  blades,  blade  twist,  and 
aspect  ratio. 

3.  The  radial  coordinates  (near  wake)  are  accurately  represented 
by  the  equation 

r  =  0.78  ♦  0.22  e 

where  A  is  a  function  of  CT  alone  and  is  given  in  Figure  62. 
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Vortex  Sheet  Axial  Coordinates 


As  discussed  in  the  subsection  entitled  Flow  Visualiation  Data,  it 
was  convenient  to  approximate  the  cross  sections  of  the  vortex  sheet  shed 
by  the  inboard  portions  of  the  blades  as  varying  linearly  with  r.  These 
linear  approximations  are  extended  (as  in  Figure  10)  at  both  ends  until 
they  intercept  the  axis  of  rotation  (r  =  0)  or  an  imaginary  cylinder  of 
radius  ?  =  1.  The  axial  coordinates  of  the  intercepts  were  found  to 
be  approximated  well  by  the  following  equations : 

f(K'r=,)*w  for  o  <  <^w  <  — 


(K2r=l)(V^r),0r  V 


2  7 r 
b 


•(4) 


o  for  o  <  <  f 

2-  .  =  i  (5) 

*r=o 

f)  ,or  ‘/'w2  -f 

The  fairings  through  the  data  shown  in  Figure  13  are  based  on  these 
equations.  Although  better  fits  to  the  data  for  r  =  0  could  be  obtained 
for  individual  conditions,  the  parameters  used  in  Equation  (5)  above 
appeared  to  give  the  best  overall  fit  to  the  r  =  0  data  for  the  conditions 
analyzed.  Further  refinements  were  not  considered  valid  in  view  of  the 
scatter  present. 


An  attempt  was  made  to  express  the  constants  appearing  in  Equations 
(4)  and  (5)  in  terms  of  more  fundamental  rotor  parameters,  as  was  done 
for  the  tip  vortex  coordinates.  Analysis  of  the  results  indicated  that 
reasonable  correlation  of  the  data  with  rotor  thrust  coefficient  existed, 
as  shown  in  Figures  63  through  65.  There  appeared  to  be  little  measurable 
effect  of  twist,  aspect  ratio,  tip  speed,  or  numbers  of  blades  on  the 
outboard  axial  coordinate  parameters,  K^-  Q  and  K2-  .  (Figures  63 

and  64).  However,  some  effect  of  twist  was  noted  on  tfie  inboard  axial 
coordinate  parameter,  Kg—  =  Q,  as  shown  in  Figure  65.  Although  scatter 
is  present,  the  data  of  Figure  65  show  a  reasonably  consistent  trend  of 
increasing  Kg~  _  Q  with  increasing  twist.  This  means  that  the  axial 
velocities  of  the  air  through  the  inner  portion  of  the  rotor  increase 
with  twist.  This  appears  to  be  reasonable  since  twist  is  used  to  increase 
the  loading  and,  hence,  the  circulation  of  the  inboard  blade  sections. 
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With  the  social  locations  of  the  inboard  sheet  cross  sections 
established  (from  Figures  63  through  65)  together  with  the  radial 
coordinates  of  the  tip  vortex  (from  Figures  58  through  62  or  Equation 
(3)),  the  radial  coordinates  of  the  inboard  vortex  filaments  can  be 
computed  using  Equation  (l).  This  equation  essentially  assumes  that  the 
contraction  of  the  inboard  vortex  filaments  (as  measured  along  the  sheet 
cross  section)  is  determined  by  the  degree  of  contraction  of  the  tip 
vortex  at  the  axial  location  where  the  inboard  sheet  extension  intersects 
the  tip  vortex  trajectory. 


Generalized  Wake  Equations 


The  generalized  wake  equations  (applicable  in  the  stable  wake  region) 
are  summarized  below. 


1.  Tip  Vortex  Axial  Coordinates: 


k|  \f/w  for  0  <  \j/w  <  k 


<*T>  V  is  +  *2  (*„  -  x)  for  **  2  nf 

k|  -  -0.25  (cT/a  +  0.001  g|  ) 

-(|.4I  +  0.0141  -JcTfc 


~-(l+00'  0'de,)^ 

2.  Tip  Vortex  Radial  Coordinates: 

r  =  A  +  ( 1  -A)e"X^w 
A  =  0.78  (neor  wake) 
X  =  0.145  +  27  CT 

3.  Vortex  Sheet  Axial  Coordinates: 


(Mf.i 


V,  K  ,or  0  5  s  — 

V,  '  - 22  TV2 

0  for  0  <  \jf  <  -j- 

!?  0  1  w  c 

(K*r.o)(*w‘  f)  for*w*  f 
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4.  Vortex  Sheet  Radial  Coordinates: 

(see  Equation  (1)) 

It  is  recognized  that  some  of  the  above  equations  and  constants  for  the 
wake  parameters  depend  on  the  manner  in  which  the  data  is  faired.  The 
fairings  used  represent  our  interpretation  of  the  data  and  are  believed 
to  be  consistent  and  reasonable,  taking  into  consideration,  as  they  do, 
all  of  the  data  available  for  the  different  blade  dpsigns.  However,  it 
will  be  shown  that  the  analytical  rotor  performance  results  are  very 
sensitive  to  wake  geometry.  In  view  of  this  sensitivity,  small  systematic 
refinements  to  the  fairings  and  thus  the  above  wake  constants  may  be 
required  if  consistently  accurate  performance  predictions  are  to  be 
achieved . 


COMPARISONS  OF  TIP  VORTEX  GEOMETRY  WITH  OTHER  SOURCES 

The  model  rotor  tip  vortex  geometry  was  compared  with  available  data 
from  other  sources  to  investigate  the  consistency  between  data  obtained 
on  different  rotors,  on  different  test  stands,  and  by  different  personnel. 
The  available  data  was  found  to  be  quite  limited  in  that  no  one  had 
previously  conducted  a  test  program  such  as  the  one  described  herein,  in 
which  wake  geometry  was  obtained  from  systematic  variations  of  the  blade 
design  and  operating  parameters. 

One  of  the  most  important  considerations  concerning  model  rotor  wake 
geometry  is  whether  it  agrees  with  that  of  full-scale  rotors.  Tip  vortex 
geometry  obtained  from  the  modal  rotor  tert  was  compared  with  limited 
full-scale  CH-53A  data  measured  on  the  Sikorsky  whirl  stand  as  presented 
by  Clark  and  Leiper  in  Reference  9.  A  comparison  of  the  model  and  full- 
scale  tip  vortex  coordinates  is  shown  in  Figure  66.  The  model  coordinates 
were  obtained  by  extrapolating  the  generalised  model  rotor  wake  to  a 
thrust  coefficient  of  0.01.  These  coordinates  are  in  good  agreement  with 
the  full-scale  tip  vortex  coordinates.  If  this  correlation  is  represen¬ 
tative,  it  may  be  concluded  that  the  generalised  wake  determined  from 
model  testing  is  applicable  to  full-scale  rotors.  However,  since  it  will 
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THEORETICAL  METHODS  FOR  PREDICTING  HOVER  PERFORMANCE 


A  principal  objective  of  this  investigation  was  to  evaluate  the 
accuracy  of  certain  theoretical  methods,  including  several  developed  at 
UARL,  for  predicting  rotor  hover  performance  —  particularly  as  affected 
by  assumptions  made  regarding  rotor  wake  geometry.  The  methods  considered 
in  this  study  are  listed  below. 

1.  Blade  Element  -  Momentum  Analysis 

2.  Golds tein-Lock  Analysis 

3.  Prescribed  Classical  Wake  Analysis 

4.  Prescribed  Experimental  Wake  Analysis 

5.  Prescribed  Theoretical  Wake  Analysis 

Each  method  is  described  below  with  emphasis  placed  on  the  last  three, 
inasmuch  as  the  first  two  methods  have  been  in  widespread  use  for  many 
years. 


BLADE  ELEMENT  -  MOMENTUM  ANALYSIS 


As  described  briefly  in  Reference  1,  this  analysis  is  based  on  the 
assumption  that  the  lift  acting  on  an  annulus  of  the  rotor  disc  is  equal 
to  the  change  in  momentum  of  the  air  passing  through  that  annulus.  Each 
annulus  or,  equivalently,  each  blade  section  is  assumed  to  operate 
independently  of  all  other  sections.  The  relations  developed  can  be  shown 
(Heyson,  Reference  15)  to  be  equivalent  to  those  obtained  using  vortex 
theory  in  which  the  rotor  is  modeled  by  an  infinite  number  of  blades  and 
the  vorticity  deposited  in  the  wake  of  the  rotor  forms  a  continuous 
cylindrical  vortex  sheet  having  a  diameter  equal  to  the  rotor  diameter. 

The  equations  relating  local  blade  thrust  and  local  induced  velocity  at 
the  disc  are  solved  iteratively  on  a  digital  computer  using  appropriate 
two-dimensional  airfoil  data  to  account  for  any  stall  or  compressibility 
effects  that  may  be  present.  Losses  due  to  flow  around  the  tips  of  the 
blades  are  accounted  for  by  specifying  a  ’’tip  loss  factor”,  which  assumes 
complete  loss  of  lift  over  a  small,  arbitrary  percentage  of  the  blade  tip 
region.  As  used  herein,  the  parameters  specified  as  input  to  the  analysis 
were  rotor  radius,  solidity,  tip  speed,  blade  twist,  collective  pitch, 
blade  coning  angle,  air  density,  speed  of  sound,  airfoil  c^  and  cd  data, 
and  tip  loss  factor  (0.97) • 
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THEORETICAL  METHODS  FOR  PREDICTING  HOVER  PERFORMANCE 


A  principal  objective  of  this  investigation  was  to  evaluate  the 
accuracy  of  certain  theoretical  methods,  including  several  developed  at 
UARL,  for  predicting  rotor  hover  performance  —  particularly  as  affected 
by  assumptions  made  regarding  rotor  wake  geometry.  Die  methods  considered 
in  this  study  are  listed  below. 

1.  Blade  Element  -  Momentum  Analysis 

2.  Goldstein-Lock  Analysis 

3.  Prescribed  Classical  Wake  Analysis 

4.  Prescribed  Experimental  Wake  Analysis 

5.  Prescribed  Theoretical  Wake  Analysis 

Each  method  is  described  below  with  emphasis  placed  on  the  last  three, 
inasmuch  as  the  first  two  methods  have  been  in  widespread  use  for  many 
years. 


BLADE  EI£M£NT  -  MOMENTUM  ANALYSIS 


As  described  briefly  in  Reference  1,  this  analysis  is  based  on  the 
assumption  that  the  lift  acting  on  an  annulus  of  the  rotor  disc  is  equal 
to  the  change  in  momentum  of  the  air  passing  through  that  annulus.  Each 
annulus  or,  equivalently,  each  blade  section  is  assumed  to  operate 
Independently  of  all  other  sections.  The  relations  developed  can  be  shown 
(Hey son,  Reference  15)  to  be  equivalent  to  those  obtained  using  vortex 
theory  in  which  the  rotor  is  modeled  by  an  infinite  number  of  blades  and 
the  vortlclty  deposited  in  the  waks  of  the  rotor  forms  a  continuous 
cylindrical  vortex  sheet  having  a  diameter  equal  to  the  rotor  diameter. 

The  equations  relating  local  blade  thrust  and  local  induced  velocity  et 
the  disc  are  solved  iteratively  on  a  digital  computer  using  appropriate 
two-dimensional  airfoil  data  to  account  for  any  stall  car  compressibility 
effects  that  may  be  present.  Losses  due  to  flow  around  the  tips  of  the 
bledes  are  accounted  for  by  specifying  e  "tip  loss  factor**,  which  assumes 
complete  loss  of  lift  over  e  smell,  arbitrary  percentage  of  the  blade  tip 
region.  As  used  herein,  the  peramsters  specif  ed  as  input  to  the  analysis 
were  rotor  radius,  solidity,  tip  speed,  blade  twist,  collective  pitch, 
blade  coning  angle,  air  density,  speed  of  sound,  eirfoi]  c{  and  c4  data, 
and  tip  loss  factor  (0.9?). 
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In  suomary,  the  Blade  Element  -  Momentum  Analysis  neglects  effects 
due  to  the  finite  number  of  blades  as  well  as  those  related  to  wake 
contraction. 


GOLDSTEIN-LOCK  ANALYSIS 


This  analysis  is  effectively  the  rotary-wing  equivalent  of  the 
classical  lifting-line  analysis  used  successfully  for  fixed  wings.  The 
analysis  is  based  principally  on  the  work  of  Goldstein  (Reference  l6), 
who  obtained  a  solution  for  the  velocity  potential  for  the  flow  about  an 
axially  translating,  doubly- infinite,  rigid  helicoidal  surface.  This 
surface  was  shown  by  Betz  (Reference  17)  to  represent  the  minimum  energy- 
wake  of  a  propeller  (or  rotor)  having  a  finite  number  of  blades.  By 
satisfying  the  flow  boundary  conditions  on  the  helical  surfaces 
representing  the  wake,  the  optimum  (or  Goldstein)  distribution  of 
circulation  in  the  wake  was  obtained.  Goldstein’s  results  were  applied  by 
Lock  (Reference  18),  who  showed  how  the  results  could  be  used  to  design 
a  propeller.  This  wa*  accomplished  by  assuming  the  propeller  to  be  lightly 
loaded  so  that  the  wake  was  essentially  uncontracted,  in  which  case  the 
circulation  distribution  in  the  wake  can  be  related  directly  to  the 
circulation  distribution  on  the  blade.  Through  this  assumption  and  the 
use  of  the  Goldstein  velocity  potential  to  define  the  induced  velocities 
at  the  plane  of  the  propeller,  the  blade  twist  and  chord  distributions 
necessary  tc  produce  the  Goldstein  (optimum)  circulation  can  be  determined. 
Lock  also  postulated  techniques  for  handling  the  inverse  problem,  wherein 
the  blade  geometry  is  defined  and  the  circulation  distribution  and  the 
associated  propeller  performance  are  required.  Ihis  situation  arises, 
for  example,  when  one  designs  an  optimum  propeller  for  operation  at  a 
specified  design  condition  and  then  wishes  to  know  the  performance  of  this 
propeller  at  off-desigu  operating  conditions.  In  his  approach.  Lock 
assumed  that  the  circulation  at  each  blade  seetion  was  part  of  a  Goldstein 
optimum  circulation  distribution  but  allowed  the  optimum  distribution 
associated  with  each  section  to  be  different.  This  implies  a  different 
wake  pitch  angle  variation  with  radius  than  would  be  the  case  for  a  blade 
whose  local  circulation  values  formed  part  of  the  same  optimum  distribu¬ 
tion.  The  resulting  wake  structure  is,  therefore,  technically  inconsis¬ 
tent  with  the  optimum  wake  assumptions  made  by  Goldstein;  however, 
reasonable  answers  are  expected  for  conditions  where  large  departures 
from  the  Goldstein  optimum  circulation  distribution  are  not  involved. 
Numerical  implementation  of  the  Goldstein-Lock  method,  which  involves  an 
iteration  at  each  blade  section  between  the  local  circulation  and  local 
wake  pitch  angle,  was  accomplished  using  an  existing  computer  program 
provided  by  the  Sikorsky  Division. 
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In  suan&ry,  the  Goldstein-Lock  Analysis  accounts  for  the  finite 
number  of  blades  on  a  rotor  (thereby  eliminating  the  need  for  fictitious 
tip  loss  factors)  but  still  retains  the  assumptions  that  the  blades  are 
lifting  lines  and  that  the  wake  is  uncontracted  (light  loading). 


UARL  PRESCRIBED  WAKE  HOVER  PERFORMANCE  PROGRAM 


The  next  three  analyses  to  be  discussed  all  employ  the  UARL 
Prescribed  Wake  Hover  Performance  Program  (or  briefly,  the  Prescribed  Wake 
Program)  for  the  solution  of  the  blade  circulation  and  inflow  distribution 
and  the  corresponding  integrated  rotor  performance.  Complete  generality 
(within  the  framework  of  the  assumptions  to  be  mentioned)  regarding  the 
specification  of  the  geometry  of  the  wake  was  maintained  in  the  computer 
program.  This  generality  permits  the  evaluation  of  a  wide  variety  of  wake 
geometries  such  as  a  classical  uncontracted  wake  geometry  (hereafter 
referred  to  as  the  classical  wake),  a  realistic  model  of  the  experimental 
wake,  and  a  theoretically  defined  wake.  The  incorporation  of  these  three 
wake  models  in  the  Prescribed  Wake  Program  results  in  the  analyses  11  , ted 
below ; 

1.  Prescribed  Classical  Wake  Analysis 

2.  Prescribed  Experimental  Wake  Analysis 

3.  Prescribed  Theoretical  Wake  Analysis 

These  analyses  are  identical  except  for  the  representation  of  the  wake. 

A  description  of  the  wake  model  used  in  each  will  be  presented  in  the 
following  sections.  Other  than  assumptions  regarding  wake  shape.  +he 
following  are  the  major  assumptions  that  currently  exist  in  the  Prescribed 
Wake  Program: 

1.  Each  blade  is  represented  by  a  lifting  line  (bound  vortex) 
divided  into  a  finite  number  of  segments  (blade  segments),  eaeh 
having  a  different  circulation  strength.  The  aerodynamic 
characteristics  determined  at  the  centers  of  each  segment  are 
assumed  to  be  representative  of  the  entire  segment. 

2.  The  wake  is  represented  by  a  finite  number  of  vortex  filaments 

trailing  from  the  blade  segment  boundaries.  Each  filament  is 
divided  into  straight  segments,  the  lengths  of  which  are 
determined  by  a  specified  wake  asimuth  interval,  The 

circulation  strength  of  each  trailing  vortex  filament  is 
constant  along  its  length  and  is  equivalent  to  the  difference 
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between  the  circulation  values  of  its  adjacent  bound  vortex 
segments  in  accordance  with  Helmholtz  laws  of  conservation  of 
vorticxty. 

3.  'Hie  blade  and  wake  characteristics  are  assumed  to  be  independent 
of  azimuth  position. 

4.  The  airflow  at  the  blades  is  assumed  to  be  two-dimensional 
(radial  induced  velocity  components  are  neglected).  For  the 
rotor  performance  calculations,  tabulated  two-dimensional  air¬ 
foil  data  (c,  ,  Cd,a)  are  provided  which  include  compressibility 
effects  (Mach  number  variations).  For  the  circulation  calcula¬ 
tions,  a  set  of  lift  curve  slopes  and  stall  angles  of  attack 
are  provided  which  vary  with  Mach  number. 

5.  Tangential  induced  velocity  components  are  neglected. 

6.  Small-angle  assumption^  are  deluded  for  the  inflow  angles  in 
the  circulation  solution. 

7.  Following  the  blade  circulation  and  inflow  solution,  convention¬ 
al  strip  theory  is  assumed  applicable  to  compute  the  rotor 
performance  characteristics. 

The  method  differs  from  that  developed  by  Rorke  and  Wells  (Reference  8) 
in  that  the  blade  inflow  distribution  is  determined  completely  by  the 
induced  effects  of  the  wake  (by  application  of  the  Biot-Savart  law)  as 
opposed  to  introducing  approximate  momentum  considerations.  A  flow 
diagram  showing  the  reouired  input,  sequence  of  major  operations,  and  out¬ 
put  of  the  program  is  presented  in  Figure  71.  As  indicated  in  this 
figure,  the  program  is  divided  into  three  independent  parts.  The  first 
transforms  the  wake  geometry  input  to  wake  coordinates.  The  second 
contains  the  computation  of  the  wake  influence  coefficients  at  the  blade, 
a.i  defined  by  the  Biot-Gavart  law,  and  the  numerical  procedures  for 
solving  the  circulation  matrix  and  associated  induced  velocity  distribu¬ 
tion.  In  the  third  part,  performance  characteristics  are  computed.  A 
provision  for  automatic  plotting  of  the  wake  filaments  is  included. 

Sample  computer  plots  are  shown  in  Figure  72  for  a  typical  experimental 
wake  and  a  classical  wake  model.  In  the  Prescribed  Theoretical  W’ke 
Analysis,  the  theoretical  wake  geooety  is  determined  from  a  separate 
program  (Wake  Geometry  Program)  and  is  used  as  input  to  the  Prescribed 
Wake  Program.  The  computer  time  required  by  the  Prescribed  Wake  Program 


normally  varies  from  approximately  15  seconds  to  2  minutes  ( UNl'VAC  1108 
computer)  depending  on  the  number  of  blades,  number  of  wake  elements,  and 
number  of  internal  iterations  required. 

For  all  applications  reported  herein,  the  blade-wake  model  in  the 
Prescribed  Wake  Program  was  represented  as  follows: 


Number  of  blade  segments  15 

Number  of  wake  filaxaents  per  blade  16 

Number  of  wake  revolutions  11 

Wake  azimuth  increment  30  deg 


The  15  blade  segments  were  distributed  such  that  10  were  spaced  at  0.02  R 
intervals  over  the  outer  20$  of  the  blade  span.  Wake  revolutions  beyond 
the  eleventh  were  found  to  have  a  negligible  effect  on  rotor  performance. 


PRESCRIBED  CLASSICAL  WAKE  ANALYSIS 


This  analysis  is,  in  many  respects,  similar 
Analysis  described  above  in  that  a  finite  number 
each  blade  is  represented  by  a  lifting  line,  and 
geometry  is  prescribed.  The  primary  differences 
follows : 


to  the  Gold stein-Lock 
of  blades  is  assumed, 
an  uncontracted  wake 
in  the  analyses  are  as 


1.  The  helical  sheets  of  vorticity  representing  the  blade  wakes  are 
approximated  by  a  finite  number  of  discrete  trailing  vortex 
filaments  to  facilitate  numerical  solution  (on  a  computer)  of 
basically  the  same  equations  which,  for  the  optimum  case, 
Goldstein  solved  analytically.  The  availability  of  a  numerical 
solution  also  permits  a  more  direct  solution  of  the  inverse 
problem  wherein  the  geometry  of  the  rotor  is  specified  as 
opposed  to  the  circulation  distribution. 

2.  The  wake  geometry  assumed  in  this  version  of  the  analysis 
differs  from  that  for  the  Goldstein  optimum  wake  in  that  the 
axial  transport  velocity  of  each  vertex  element  in  the  wake  is 
constant  with  radius  and  is  equal  to  the  momentum  value. 
Tangential  transport  velocities  are  assumed  to  be  zero.  For  the 
low  helical  wake  pitch  angles  associated  with  helicopter  rotor 
disc  loadings,  the  outer  portion  of  the  wake  used  approximates 
the  Goldstein  wake. 


A  sample  plot  of  the  classical  wake  trajectory  was  presented  in  Figur-  7?. 
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The  operating  conditions  for  which  rotor  perforraance  was  computed 
ing  the  Prescribed  Classical  Wake  Analysis  are  listed  in  Table  III. 


TABLE  Til. 

PRESCRIBED  CLASSICAL  WAKE  CONDITIONS 

Rotor 

Twist 

(deg) 

Aspect 

Ratio 

No.  of 

Blades 

Tip  Speed 
(fps) 

Collective 

Pitch 

(deg) 

Model 

-.3 

18.2 

2 

700 

6,8,10 

Model 

-8 

13.2 

4 

700 

6,8,10 

Model 

-8 

18.2 

6 

700 

6,8,10 

Model 

-a 

18.2 

8 

700 

8,10 

Mode  1 

-8 

13.2 

2 

525 

6,3,10 

Model 

-8 

18.2 

6 

525 

6,8,10 

Model 

0 

18.2 

0 

700 

6,8,10 

Model 

0 

18.2 

4 

70C 

6,8,10 

Model 

0 

18.2 

6 

700 

6,8,10,11 

Model 

0 

18.2 

8 

700 

8,10 

Model 

-lo 

l9.2 

2 

700 

6,8,10 

Motel 

-16 

18. ' 

6 

700 

6,8,10 

Model 

-a 

13.'' 

2 

700 

6,8,10 

Model 

13.'-' 

o 

700 

6,8,10 

•  ■  _  >  » 

-i 1 

17. « 

•■> 

713 

7,o,ll 

C  U 

-r. 

1‘  .7 

i  < 

<^o8 

3.10,17 

A I !  -  Tid 

d*.  Ur.,;  inve 

i  gad.ed 

for1'*.;  pool 

”.o  oper*ti 

on . 

V.- 


PRESCRIBED  EXPERIMENTAL  WAKE  ANALYSIS 


This  analysis  differs  from  the  previous  analyses  in  that  a  more 
realistic,  contracted  wake  geometry  based  on  experimental  flow  visualiza¬ 
tion  data  is  used  as  input  to  the  Prescribed  Wake  Program.  Rather  than 
input  coordinates  for  each  vortex  element  in  the  wake,  Equations  (1) 
through  (5)  are  used  to  define  the  wake  and  the  individual  coordinates 
are  computed  in  the  program.  In  this  manner  it  is  only  necessary  to  input 
the  following  wake  constants  and  certain  wake  azimuth  angles  : 


A,  A,  k^,  k2,  K1?  _  0*  _  o’  Klf  =  1’  *2j  -  i 


The  required  wake  azimuth  angles  are  simply  the  junction  angles  that 
bound  the  kp  and  k2  regions.  Use  of  the  wake  equations  and  the  above 
constants  greatly  simplifies  the  input  requirements  while  retaining 
sufficient,  accuracy  for  the  computation  of  induced  velocities  at  the 
rotor  blades. 

tigure  72  indicates  how  the  computer  transform;  the  input  wake 
constants  to  coordinate  form  and  plots  the  resulting  wake  pattern.  On 
the  1,  ft  side  of  this  figure  are  the  top  and  side  views  of  the  computer 
representation  of  a  typical  experimental  wake.  For  clarity,  only  the 
wakes  from  1  blade  and  12  vortex  filaments  are  shown.  For  this  test 
condition,  5  cf  the  16  vortex  filaments  over  the  outer  8^  of  the  blade 
were  used  to  represent  the  tip  vortex.  The  spanwise  division  between  the 
vortex  sheet  and  tip  vortex  regiors  is  determined  by  the  requirement  that 
the  vortex  filament;  grouped  in  the  tip  vortex  have  the  same  circulation 
sense  and  one  which  is  consistent  with  a  negative  derivative  of  the  final 
computed  bound  circulation  distribution  (-df/dr)  over  the  tip  region  of 
the  blade.  This  results  in  the  radial  location  of  the  peak  bound 
circulation  as  the  dividing  point  between  the  inboard  sheet  and  tip  vortex 
portions  of  the  wake.  An  iteration  is  built  into  the  computer  program  to 
insure  this  consistency.  A  program  refinement,  included  under  this 
contract,  was  a  ,  improvement  of  the  wake  geometry  representation  through 
incorporation  of  a  provision  for  approximating  the  roll-up  of  the  tip 
vortex  filaments  into  a  single  filament.  This  was  accomplished  by 
truncating  the  inner  tip  vortex  filaments  at  an  input  azimuth  (30  degrees 
was  used  in  this  investigation)  and  assigning  tne  experimental  tip  vortex 
geometry  and  the  combined  circulation  strength  to  the  remaining  filament 
(Figure  72).  In  addition,  an  improved  representation  of  the  vortex  sheet 
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war,  included  which  allows  the  vortex  sheet  to  extend  to  the  tip  vortex 
boundaiy.  Finally,  a  provision  was  included  for  irputting  tip  vortex 
coordinates  rather  than  using  the  curve  fit  equations,  if  such  a  mode  of 
operation  is  desired.  The  far  wake  model  was  assumed  to  be  a  smoothly 
contracting  extension  of  the  near  wake  as  shown  in  Figure  72.  The  inclu¬ 
sion  of  provisions  accounting  for  the  instability  of  the  fir  wake  observed 
in  this  investigation  was  beyond  the  scope  of  this  3tudy.  The  exact 
nature  of  this  instability,  and  whether  the  unstable  region  is  sufficiently 
removed  from  the  rotor  so  as  not  to  affect  rotor  performance  remains  to  be 
investigated , 

The  operating  conditions  for  which  rotor  performance  was  computed 
using  the  Prescribed  Experimental  Wake  Analysis  are  listed  below. 


TABLE  IV. 

PRESCRIBED  EXPERIMENTAL  WAKE  CONDITIONS 

Rotor 

Twist 

(deg) 

Aspect 

Racio 

No.  of 

Blades 

Tip  Speed 
(fps) 

Collective 

Pitch 

(deg) 

Model 

-8 

18.2 

2 

700 

6,8,10 

Model 

-8 

18.2 

6 

700 

6,8,10 

Model 

-8 

18.2 

2 

525 

6,8,10 

Model 

-8 

18.2 

6 

525 

6,8,10 

Model 

0 

18.2 

2 

700 

8,10 

Model 

0 

18.2 

6 

700 

8,10 

Model 

-8 

13.6 

2 

7C0 

3,10 

HU-1A 

-12 

17.3 

2 

713 

7.9,11 

CH-53-A 

— -  _ 

-6 

16.7 

6 

698 

8,10,12 

All  conditions  investigated  correspond  to  CGE  operation. 


PRESCRIBED  THEORETICAL  WAKE  ANALYSIS 


General  Approach 

The  final  performance  method  evaluated  was  the  so-called  Prescribed 
Theoretical  Wake  Analysis.  In  this  method  the  wake  is  again  prescribed  as 
input  to  the  Prescribed  WcVe  Program  in  order  to  compute  the  blade  and  wake 
circulation  distributions  and  the  associated  integrated  rotor  performance. 
In  this  instance,  however,  the  input  wake  geometry  is  determined  princi¬ 
pally  from  theoretical  considerations  rather  than  from  experiment. 
Conceptually,  the  computation  of  the  theoretical  wake  geometry  can  be 
accomplished  by  the  following  steps  : 

1.  Estimate  the  vorticity  (circulation  strength)  in  the  wake  from 
a  previous  solution  of  the  bound  circulation  distribution  on 
the  blade. 

2.  Specify  an  initial  wake  geometry. 

3.  Apply  the  classical  Biot-Savart  relation  to  compute  the 
velocities  induced  in  the  wake  by  the  wake  vorticity. 

4.  Integrate  these  velocities  over  a  small  increment  in  time  to 
define  a  new  wake  geometry. 

5.  Repeat  Steps  3  and  4,  alternately,  until  a  converge!  wake 
geometry  corresponding  to  the  initial  estimate  of  blade  bound 
circulation  is  obtained. 

6.  Compute  a  new  estimate  of  the  blade  bound  circulation 
distribution  using  the  Prescribed  Wake  Program  and  the  wake 
geometry  from  Step  5. 

7.  Repeat  Steps  2  through  6,  iterating  until  a  compatible  geometry- 
circulation  solution  is  obtained. 

The  process  just  described  is  illustrated  schematically  in  figure  73-  In 
the  top  half  of  this  figure,  the  fundamental  iteration  between  the  blade 
circulation  and  wake  geometry  is  indicated,  while  the  lower  half  of  the 
figure  indicates  the  second  iteration  required  to  obtain  a  geometry 
consistent  with  the  current  estimate  of  the  blade  circulations.  A  computer 
program  termed  the  Wake  Geometry  Program  has  been  developed  to  per f era 
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tnis  second  iteration.  Iterations  of  the  type  described  are  necessary 
because  the  complexity  of  the  rotor  wake  geometry  precludes  a  closed-form 
solution  to  the  problem. 

Although  the  process  described  above  is  conceptually  straightforward, 
its  actual  numerical  implementation  represents  a  formidable  task,  even 
with  modern  computers,  because  of  the  great  number  of  vortex  elements  that 
conceivably  could  be  used  to  represent  the  wake.  Approximations  are, 
therefore,  necessary  to  limit  computing  time  requirements.  The  following 
sections  describe  the  procedures  and  approximations  employed  in  this 
study. 

Before  proceeding  to  th’.  detailed  discussion  of  the  procedures  and 
approximations  used,  a  few  general  remarks  are  in  order  regarding  the 
computer  program  developed.  The  core  program  available  for  this  study 
was  basically  one  which  was  capable  of  computing  wake  geometries  for 
forward  flight  conditions  (see  Reference  2).  As  a  result,  the  program  is 
not  optimized  for  the  hover  condition  ar.d  thus,  for  example,  the  symmetry 
features  of  the  hover  flight  condition  are  not  included.  Incorporation  of 
such  features  could  significantly  reduce  computing  time.  Also,  initial 
attempts  to  apply  the  program  in  hover  disclosed  an  apparent  instability 
of  the  wake  at  moderate  distances  from  the  rotor.  This  instability 
appeared  to  conflict  with  the  more  or  less  classical  conception  of  the 
hovering  rotor  wake  as  being  a  stable,  smoothly  contracting  wake. 
Considerable  effort  was  expended  in  investigating  the  possibility  that  the 
apparent  instability  was  a  result  of  the  numerical  procedures  used;  in 
the  process,  several  modes  of  operation  of  the  program  were  developed 
before  a  conclusion  was  reached  that  the  physical  wake  is  apparently 
unstable.  The  different  operating  modes  allow  the  user  to  check  various 
assumptions  made  in  the  computation  of  the  wake  geometry.  The  resulting 
computer  program  should,  therefore,  be  considered  a  research  program 
rather  than  one  which  is  optimized  for  maximum  efficiency  in  production 
use . 

Numerical  Procedures  and  Approximations 


The  procedures  ant  approximations  employed  in  developing  a  numerical 
met hoi  for  computing  rotor  wake  geometries  for  hovering  flight  are 
described  below.  The  reader  is  also  referred  To  Reference  7  for  more 
backer  ml  material. 


Discrete  Vortex  Representation 


The  continuous  helical  sheets  of  vorticitj  trailed  by  the  various 
bisdes  of  the  rotor  are  represented  by  a  finite  number  of  discrete  trail¬ 
ing  vortex  filaments.  This  approach  is  fairly  standard,  being  used  by 
Westwater  in  Reference  19,  for  example,  to  compute  the  geometry  of  fixed- 
wing  wakes.  Figure  7k  shows  a  schematic  representation  of  a  continuous 
sheet  wake  from  one  blade  and  the  approx inat ion  to  this  sheet  made  using 
discrete  filaments.  Because  of  the  nature  of  rotor  blade  bound  circula¬ 
tion  distributions,  the  discrete  vortices  naturally  group  into  two  parts: 
a  strong,  rolled-up  tip  vortex  filament  ar.d  several  weaker  filaments 
representing  the  inboard  portion  of  the  vortex  sheet.  In  addition,  it  is 
evident  from  Figure  7k  that  the  discrete  vertices  are  each  represented  by 
a  series  of  straight-line  segments  rather  than  by  continuous  curves.  Fhe 
length  of  the  segments  is  determined  by  the  time  increment  used  in  the  in¬ 
tegration  of  the  distorting  velocities  acting  on  the  filaments  and  the 
magnitude  cf  4  he  distorting  velocities  themselves.  Each  discrete  vortex 
was  assumed  to  have  a  finite  core  radius  equal  to  0.005  R.  Velocities 
within  the  core  were  arbitrarily  set  to  zero.  More  elaborate  representa¬ 
tion  of  the  velocities  within  the  core  are  available  (e.g.,  Reference  20) 
but  were  not  used  herein  in  view  of  other,  more  important  assumptions 
made  in  this  analysis. 

Far  Wake  Representation 

Theoretically,  the  helical  wake  under  the  rotor  extends  downstream 
to  infinity.  As  indicated  previously,  the  basic  computer  program  available 
for  this  study  was  a  forward- flight  program.  In  this  program  the  wake  is 
truncated  after  a  specified  number  of  revolutions  as  indicate!  in  Figure 
75  (for  a  hover  condition).  Initial  attempts  to  compute  the  distortions 
of  such  a  truncated  wake  in  hover  showed  that  it  quickly  tended  to  roll  up 
as  shown  in  Figure  75-  In  Figure  75(a),  the  truncated  wake  at  the  start 
of  the  computation  is  shown;  while  : n  Figure  75(b),  the  early  stages  of 
roll-up  are  evident  after  12  time  steps.  This  phenomenon  is  similar  to 
that  observed  in  Reference  ?1,  where  distortions  of  a  wake  represented  by 
a  finite  number  of  vortex  rings  were  computed.  The  roll-up  observe!  is 
reasonable  inasmuch  as  symmetry  const terat ions  dictate  that  the  self- 
induced  radial  velocities  acting  on  the  truncated  helical  v  rtex  must  be 
ant i symmetric  about  the  midpoint  of  '  .e  helix.  In  an  attempt  to  av  i \ 
this  problem  (which  was  con si lere  '  be  a  fictitious  problem  cause!  -nly 
by  the  wake  truncation  feature )  an!  at  the  same  time  keep  computing  times 
low,  a  modification  to  the  basic  program  was  :ua!e  wherein  the  overall  wake 
was  divided  into  near-  an!  far -wake  eg ions  as  shown  in  Figure  7 '  .  The 


near  wake  in  allowed  to  distort  freely  under  the  influence  of  the  wake- 
induced  velocities.  The  far  wake,  however,  was  r.rtificaily  constrained  so 
that  its  diaraetei  Dp,  was  at  each  tire  step  set  equal  to  twice  the  radius 
of  the  last  vortex  element  in  the  immediately  preceding  near  wake.  The 
axial  spacing  was  set  equal  to  the  last  revolution  of  the  initial  wake 
(which  for  these  calculations  was  based  on  experimental  data).  As  will  be 
discussed  belcw,  the  partially  constrained  far  wake  moves  away  from  the 
rotor  as  the  computation  proceeds,  so  that  its  influence  continuously 
diminishes . 

Numerical  Integration 

The  numerical  computation  of  the  distorted  wake  geometry  involves 
two  basic  steps:  (1)  the  use  of  the  Biot-Savart  Law  to  compute  the  wake 
distortion  velocities  (vx  Vy^,  v^p)  produced  by  a  given  wake  geometry  and 
vorticity  at  any  point  (P)  joining  the  straight-line  vortex  elements  (see 
Figure  77),  and  (2)  the  integration  of  these  velocities  over  a  small  time 
increment  (or  tire  step)  to  establish  a  new  wake  geometry.  These  basic 
steps  are  successively  repeated  until  there  is  no  change  in  the  computed 
wake  geometry.  The  expressions  used  to  compute  and  integrate  these 
velocities  are  given  in  Appendix  I.  As  the  wake  distorts  during  the  time 
increment,  the  blades  are  allowed  to  rotate  to  new  positions.  This 
rotation  leads  to  the  generation  of  new  near-wake  elements.  As  a  result, 
the  number  of  near-wake  elerents  continually  builds  up  while  the  partially 
constrained  far  wake  moves  axially  away  from  the  rotor.  The  process  will 
be  illustrated  by  a  specific  example  in  a  late’  section.  Other  modes  of 
operation  are  possible  with  the  program  but  were  not  evaluated  in  this 
contract.  These  will  be  described  in  the  report  documenting  the  program 
itself. 

Further  Inboard  Wake  Approximations 

As  noted  above,  the  inboard  sheet  has  been  approximated  by  a  series 
of  discrete,  segmented  vortex  filaments.  If  one  were  to  consider 
computing  the  contribution  of  each  vortex  element  in  the  wake  to  the 
velocities  acting  on  all  other  elements  in  the  wake,  one  would  find  that 
the  computing  time  quickly  gets  out  of  hand.  Significant  reductions  in 
computing  time  can  be  accomplished  by  (l)  eliminating  elements  in  the 
wake,  (?.)  avoiding  the  computation  of  the  distorting  velocities  induced  by 
certain  elements  at  each  and  every  time  step,  or  (3.)  avoiding  altogether 
the  computation  of  the  distorting  velocities  induced  by  certain  elements 
or.  othe-  elements.  A  combination  of  all  three  approaches  was  used  in  this 
study  (although  it  should  be  noted  that  the  program  has  optional  modes  of 


operation  which  will  allow  a  complete  interaction  computation  if  one  so 
desires).  Thus,  the  far  wake  was  truncated  after  Np  revolutions.  Also, 
advantage  was  taken  of  the  fact  that  the  circulation  strength  of  each 
inboard  vortex  filament  is  significantly  less  than  that  of  the  tip  vort«*x. 
The  distorting  velocities  induced  by  the  inboard  elements  on  the  tip 
vortex  are,  therefore,  generally  small  compared  to  those  induced  by  the 
tip  vortex  itself.  Rather  than  neglect  the  inboard  wake  completely  in 
the  wake  geometry  computation  a.-  was  done  in  Reference  2,  it  was  decided 
to  (1)  specify  a  representative  inboard  wake-tip  vortex  geometry, 

(2)  compute  the  velocities  Induced  along  the  tip  vortex  by  the  inboard 
wake  for  this  geometry  at  the  initial  time  step,  and  (3)  keep  these 
velocities  constant  for  all  succeeding  time  steps  during  which  the  tip 
vortex  distorts  under  both  its  <wn  variable  influence,  which  varies  with 
time,  and  that  of  the  inboard  wake,  which  is  constant  with  time .  The 
representative  wake  geometry  used  to  compute  the  inboard  wake  effect  on 
the  tip  vortex  was  estimated  for  each  operating  condition  from  the 
experimental  wake  data  given  in  Figures  53  through  65.  A  typical  geometry 
is  shown  in  Figure  j8.  As  shown  in  this  figure,  the  outer  boundary  of 
the  inboard  wake  was  terminated  at  a  distance  about  0.04  R  from  the  tip 
vortex  trajectory,  except  in  the  immediate  vicinity  of  the  blade,  where, 
of  course,  the  radial  location  was  determined  from  the  nature  of  the 
bound  circulation  distribution.  The  resulting  boundary  appeared  reason¬ 
able  when  compared  with  smoke  picture  observations  and  considering  that 
the  local  trajectory  of  a  v'-  tex  element  should  be  dominated  by  the 
rotational  velocity  induced  by  the  nearest  element  of  the  tip  vortex. 
Finally,  specification  of  the  inboard  sheet  geometry  eliminated  the 
necessity  of  computing  the  effect  cf  the  tip  vortex  on  the  inboard  sheet. 
By  incorporating  these  approximations,  the  computing  times  required  were 
reduced  by  about  an  order  of  magnitude  to  about  ^  minute  per  time  step  f'Y 
a  typical  tvc-bladed  rotor  case  and  2^  minutes  per  time  step  for  a  typical 
six-bladed  rotor  ease  (UKTVAC  1106  Computer). 

Bound  Vortex  Effect 


The  velocities  induced  on  the  tip  vortex  by  the  slide  bound  vortieity 
were  neglected.  This  eliminated  an  artificial  numerical  problem  resulting 
from  the  representation  of  the  blade  as  a  lifting  line  and  from  the 
computation  of  velocities  at  discrete  points  in  ihe  wake.  As  a  result, 
unrealistically  large  fluctuations  in  induced  velocity  occurred  as  the 
bound  vortex  passed  over  points  on  the  ti>  vortex  generated  by  the  pre¬ 
vious  blade.  Removal  of  these  basic  restrictions  to  the  program  was 
beyond  the  scope  of  this  study.  Because  the  generally  antisymmetric 
nature  of  the  velocities  induced  by  the  bound  vortex  and  their  rapid  decay 
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with  increasing  distance  from  the  bound  vortex,  it  is  believed  that  the 
neglect  of  these  velocities  should  have  little  effect  on  the  computed 
average  position  of  the  tip  vortex. 

Operating  Conditions  and  Program  Parameters  Selected 

The  operating  conditions  for  which  rotor  wake  geometries  and 
associated  performance  characteristics  were  computed  using  the  Prescribed 
Theoretical  Wake  Analysis  are  given  below. 


TABLE  V. 

PRESCRIBED 

THEORETICAL  WAKE  CONDITIONS 

Rotor 

Aspect 

Ratio 

Twist, 

(deg) 

No.  of 

Blades 

Tip  Speed 
(fps) 

Collective 

Pitch 

(deg) 

Model 

10.2 

-8 

2 

700 

6,8,10 

Model 

18.2 

-8 

2 

525 

8,10 

Model 

18.2 

0 

2 

700 

8,10 

Model 

13.6 

-8 

2 

700 

8,10 

MM  el 

18.2 

-8 

o 

700 

8 

HU- LA 

17.3 

-12 

2 

713 

9 

CH-53A 

16.7 

-6 

6 

698 

10 

All  conditions  investigated  correspond  to  OGE  operation. 


In  computing  the  wake  geometries  ar.d  rotor  performance  for  the 
conditions  described  above,  the  following  parameters  were  selected  ar 
input  to  the  Wake  Geometry  Program  (input  to  the  Prescribed  Wake  Program 
was  the  .-.aoe  as  that  described  in  the  section  entitled  "UARL  PRESCRIBED 
WAKE  H'VER  PERFORMANCE;  PROGRAM" ) . 


1. 

Number  of  inboard  wake  filaments 

4  to  5  (as  required) 

2. 

Initial  number  of  near  wake 
revolutions  (tip  vortex) 

3/blade  for  b  =  2 
? /blade  for  b  =  6 

3. 

Far-wake  revolutions  (tip  vortex) 

3/blade  for  b  =  2, 

3/blade  for  b  =  6 

4. 

Number  of  time  steps  used  in 
integration 

49  for  b  =  2, 

29  for  b  *  6, 

5. 

Nondimensional  time  interval 
used  in  integration,  Ai/f=OAt 

0.524  (30  deg)  for  b  =  2, 
0.349  (20  deg)  for  b  -  6 

6. 

Number  of  inboard  wake 

revolutions 

4  for  b  =  2, 

3  for  b  =  6 

7. 

Initial  wake  geometry 

estimated  from 
figures  53  through  65. 

Although 
for  the 

the  number  of  filaments  per  blade  i 
six-bladed  rotor  because  of  program 

n  tne  initial  wake  was  less 
storage  limitations,  the  total 

number  of  filaments  in  the  wake  for  the  two-  and  six-bladed  rotors  were 
comparable. 

Wake  Stability  Characteristics 

Evidence  of  Wake  Instability 

Initial  results  for  the  two-bladed  model  rotor  obtained  from  the 
Wake  Geometry  Program  (as  modified  for  this  study)  disclosed  an  aoparent 
instability  of  the  helical  tip  vortex.  This  instability,  which  occurred 
at  axial  distances  from  the  rotor  plane  as  low  as  0.15  R,  is  illustrated 
in  Figure  79*  Here  the  computed  axial  and  radial  coordinates  of  the  tip 
vortex  for  several  different  time  steps  in  the  calculation  are  presented. 
At  the  start  of  the  computation  (extreme  left-hand  panel  of  Figure  79), 
the  wake  is  assumed  to  smoothly  contract  and  to  be  composed  of  the  near- 
and  far-wake  regions  described  previously.  That  portion  cf  the  near  wake 
in  the  immediate  vicinity  of  the  rotor  plane  was  obtained  from  t-he 
generalized  experimental  results  of  Figures  53  through  o5 .  The  remainder 
of  the  initial  wake  was  obtained  through  extrapolation  in  a  manner 
consistent  with  the  classical  conception  of  a  smoothly  contracting  rotor 
wake  in  hover.  As  the  calculation  of  the  geometry  proceeds  (i.e.,  as  the 


number  of  time  steps  increases),  the  rotor  rotates,  new  near-wake  vortex 
elements  are  created,  the  far  wake  moves  downstream,  and  the  near  wake 
distorts  under  both  its  own  influence  and  that  of  the  far  wake.  The 
remaining  panels  of  Figure  79  illustrate  the  process  by  showing  the  wake 
as  computed  after  6,  12,  lB,  and  2b  time  steps.  (For  the  nondimens ional 
time  used  in  this  computation,  0.52^,  12  time  steps  correspond  to  one 
revolution  of  the  rotor.)  The  circles  shown  on  the  wake  boundary 
represent  the  intercepts  of  the  helical  tip  vortices  of  both  blades  with 
the  r  -  z  plane  containing  one  of  the  blades  as  a  reference.  These 
circles  would  mark  the  positions  of  the  vortices  that  would  be  observed 
in  photographs  taken  in  a  flow  visualization  experiment  (such  as  conducted 
under  this  contract)  where  smoke  was  introduced  into  the  rotor  flow  field 
in  the  nonrotating  system.  It  is  evident  from  Figure  79  that  the  wake 
becomes  unstable  at  some  distance  belcw  the  rotor. 

In  view  of  the  importance  of  the  conclusion  regarding  the  stability 
of  the  rotor  wake,  additional  checks  on  the  numerical  procedures  being 
used  were  made.  In  addition,  available  literature  and  data  ve.--e  reviewed 
in  an  att  ’capt  to  obtain  supporting  evidence.  The  numerical  checks 
included  the  use  of  different  nondimensional  time  interval,  and  other 
numerical  integration  techniques  such  as  the  Runge-Kutta  method  All 
results  obtained  were  qualitatively  the  sat*'. 

The  review  of  the  literature  disclosed  an  analysis  by  Levy  and 
Forsdykc  (Reference  22)  which  indicated  that  a  aoubly-infinite ,  constant- 
diameter,  helical  vortex  can  be  unstable  when  the  pitch  angle  of  the 
helix  is  less  than  0.3  radian.  Helix  pi.ch  angles  for  helicopter  hover 
conditions  lie  in  this  critical  cangt.  Although  the  tetual  rotor  wake  in 
hover  is  a  contracting  wake,  it  is  believed  that  the  results  from 
Reference  22  imply  the  possibility  of  a  potential  stability  problem. 
Comparison  of  the  unstable  wake  results  obtained  herein  with  results  of 
other  analyses  In  which  atteu.pt  3  '-ere  mide  to  ec.aj.ate  the  geometry  of 
rotor  wakes  in  hover  was  generally  inconclusive  because  of  differing 
assumptions  and  numerical  techniques  employed.  For  example,  interpre¬ 
tation  of  the  results  of  Reference  21  (where  vortex  rings  are  used  and  the 
equations  are  integrated  as  in  this  studv)  are  complicated  by  the  general 
roll-up  problem  of  a  tr mcated  wake  noted  earlier.  However,  a  lack  of 
convergence  of  certain  areas  of  tiie  wake  was  noted.  In  Reference  9,  an 
analysis  employing  helical  vortex  filaments  but  a  fundame  tally  different 
numerical  computation  scheme  is  described.  Also,  only  two  fully  distor-- 
ing  revolutions  of  the  wake  are  considered  together  with  a  pari .ally  eon- 
stiained  fa*  wake.  Qualitative  convergence  tests  ure  employed,  and  no 
mention  of  any  observe  1  wake  instability  is  made  in  connection  with  the 


limited  results  presented.  In  spite  of  apparent  differences  regarding 
wake  stability,  results  of  Reference  9  appear  to  agree  well  with  results 
obtained  in  this  study  for  the  region  of  the  wake  near  the  rotor  plane. 

This  will  be  discussed  further  in  a  subsequent  section. 

Experimental  flow  visualization  data  reviewed  included  that  from 
References  1,  10,  and  12,  as  well  as  data  obtained  under  this  contract. 

The  review  led  to  the  following  observation  : 

1.  The  critical  area  of  the  wake  where  the  instability  was 
predicted  was  generally  poorly  defined. 

2.  In  a  large  number  of  cases,  evidence  of  a  rolling  of 
successive  coils  of  the  helical  tip  vortexes  around  one 
another  (local  roll-up)  could  be  found.  In  those  experiments 
where  smoke  was  injected  into  the  nonrotating  system,  such  as 
the  present  study,  this  roll-up  was  inferred  from  an  uneven 
axial  spacing  of  the  tip  vortexes  and,  more  importantly,  from 
the  fact  that  downstream  vortexes  had  radial  coordinates 
greater  than  those  for  upstream  vortexes. 

3.  Shortly  after  evidence  of  roll-up  of  the  experimental  wake  was 
noted,  further  tracking  or  observation  of  the  tip  vortexes 
became  exceedingly  difficult.  Thus,  the  roll-up  may  be  the 
cause  of  the  poorly  defined  wake  noted  in  Item  1  above. 

k.  In  no  case  was  a  smoothly  contracting  tip  vortex  observed  over 
large  enough  axial  distances  below  the  rotor  as  to  definitely 
preclude  the  possibility  of  an  instability. 

These  observations  are  substantiated  by  the  sample  photographs 
presented  previously  (Figure  30)  and  the  additional  photographs  presented 
in  Figures  80  through  82.  A  schematic  interpretation  of  these  photographs 
showing  local  roll-up  of  the  tip  vortex  is  shown  in  Figure  83.  In  summary, 
then,  the  prediction  of  an  unstable  wake  does  not  appear  to  be  inconsis¬ 
tent  with  available  evidence;  further  work  should,  however,  be  undertaken 
to  define  the  characteristics  of  experimental  wakes  in  more  detail  in  the 
critical  areas.  'tandard  flow  visualization  techniques  appear  to  be  of 
marginal  usefulness  of  this  purpose,  and  hot-wire  and/or  holographic 
approaches  should  be  considered. 
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Inpli.c3ti.ons  of  Instability  on  Procedures  Employed 


The  existence  of  a  wake  instability  obviously  complicates  both  the 
computation  of  wake  geometries  and  the  use  of  these  geometries  in  the 
subsequent  computation  of  rotor  performance.  Some  modification  of 
anticipated  procedures  was  necessitated  and  those  procedures  that  were 
used  are  described  in  the  following  paragraphs. 

Close  examination  of  the  geometry  results  obtained  indicated  that 
although  the  computed  wake  did  become  unstable  at  moderate  distances  from 
the  rotor,  the  portion  of  the  wake  in  the  immediate  vicinity  of  the  rotor 
plane  (extreme  near-wake)  did  converge  to  a  reasonably  stable  answer, 
with  the  degree  of  convergence  improving  with  increasing  proximity  to  the 
rotor.  This  presumably  is  due  to  the  fact  that  one  end  of  the  wake  is 
perfectly  stable,  being  tied  as  it  is  to  the  blade  itself.  The  conver¬ 
gence  of  the  extreme  near-wake,  which  was  also  observed  in  the  study  of 
Reference  21,  is  illustrated  in  the  sample  results  presented  in  Figures 
R-*  and  8b.  In  Figure  34,  the  time  histories  of  the  radial  and  axial 
coordinates  for  two  critical  points  in  the  wake  of  a  two-bladed  rotor  are 
shown.  These  are  the  points  of  the  wake  lying  immediately  below  the  blade 
(*v  3  l80  deg  and  3  60  deg).  In  Figure  85,  on  the  other  hand,  the  computed 
radial  and  axial  coordinates  at  two  different  time  steps  are  presented  as 
a  function  of  wake  azimuth  position.  Experimental  results  are  also  shown 
in  Figure  85  for  reference  purposes.  Convergence  of  the  extreme  near-wake 
to  a  stable  solution  is  apparent.  Also,  results  such  as  these  led  to  the 
selection  of  49  time  steps  in  the  wake  geometry  integration  as  being 
adequate  for  all  two-bladed  rotor  computations.  This  corresponds  to  a 
computation  of  the  wake  for  a  time  period  corresponding  to  four  rotor 
revolutions  or  eight  blade  passage  cycles.  For  the  six-bladed  rotor,  only 
?9  time  steps  (i.e.,  1.6l  rotor  revolutions,  9* 06  blade  passage  cycles) 
were  used  because  of  machine  storage  limitations.  It  was  also  necessary 
to  re  luce  the  time  interval  to  0.349  (?C  deg)  to  avoid  a  mathematical 
instability.  In  spite  of  the  lower  number  of  tine  steps,  reasonable  near¬ 
wake  convergence  was  obtained. 


'n  using  the  predicted  tip  vertex  ge 
circulation  Ustributiens  by  means  f  the 
c-'nver.el  near-wake  results  were  in  most 
same  analytic  expressions  :  ?  «ua“  *  ns 
<r,vr  a  1  wake  la*  a;  namely. 


one tries  to  compute  new  blade 
Prescribe!  Wake  Program,  the 
instances  approximated  by  the 
an*.  ise!  in  analyzing  the 
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Constants  in  the  expressions  were  determined  by  emphasizing  correct 
positioning  of  the  vortices  in  the  immediate  vicinity  of  the  blade  (i.e., 
atV\r=  180  deg  and  3$0  deg  for  a  two-bladed  rotor).  In  an  isolated  case, 
where  the  above  analytical  expressions  did  not  appear  to  represent  a  good 
fit  to  the  computed  geometry,  the  actual  computed  r,  z  values  of  the  extreme 
near-wake  were  used  and  extrapolated  in  a  reasonable  manner  to  obtain  the 
entire  wake.  By  using  these  procedures,  the  necessity  for  employing 
computed  wake  coordinates  for  that  portion  of  the  wake  which  was  predicted 
to  bt  unstable  was  avoided. 

Sample  Wake  Geometry-Blade  Circulati  in  Computation 

In  view  of  the  relative  complexity  of  the  procedures  involved  in 
computing  the  wake  geometry  and  associated  rotor  performance,  a  specific 
example  illustrating  the  various  steps  of  the  process  is  desirable  and  is 
presented  below. 

Rotor  and  Rotor  Operating  Condition 


Number  of  Blades  2 

Twist  -8  deg 

Tip  Speed  TOO  fps 

Aspect  Ratio  13 .0 

Nominal  Cj/tr  O.C-95 

Radius  2.2}  ft 

Collective  Pitch  10  deg 
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:T  ice  lure.: 


tmate  the  wake  gee  •  - ■;  ruing  the  experimental  results  of 

Figures  59  through  ’>5  (nee  Figure  %).* 

O-mpute  the  bound  circulation  distribution  using  the  experimental 
wake  geometry  from  .‘tep  1  as  input  to  the  Prescribed  Wake  Program, 
ihese  results  are  shown  in  Figure  *'7. 

Approximate  the  continuous  bound  circulation  distribution  from  Step  2 
by  five  constant -circulation  segments  as  shown  in  Figure  87-  Note 
that  the  circulation  peak  (  1' max )  near  the  tip  of  the  blade  has  been 
increased  in  magnitude  and  decreased  in  radial  extent  in  making  the 
approximation.  This  is  done  in  anticipation  of  the  expected  final 
circulation  distribution.  The  particular  changes  made  were  prompted 
by  the  observation  that  the  theory  generally  predicted  the  vortex 
from  the  previous  blade  to  lie  closer  to  the  rotor  and  farther 
outbonr  1  than  was  observed  experimentally.  The  step  changes  in 
circulation,  of  course,  imply  the  generation  of  trailing  vortex 
filaments  having  circulation  strengths  equal  to  the  changes  in  bound 
circulation.  The  tip  vortex  is  trailed  from  the  r  =  1  location, 
while  the  vortexes  representing  the  inboard  sheet  trail  from  f  *  0.15, 
0.75,  3.1*,  0. 75,  and  0.385  in  this  example. 

Vising  the  approximate  bound  circulation  distribution  from  Step  3  to 
specify  the  circulation  strengths  of  the  inboard  and  tip  vortex 
filaments  and  using  the  geometry  shown  in  Figure  87  for  the  initial 
starting  wake,  compute  the  corresponding  distortions  of  the  tip 
vertex  using  the  Wake  geometry  Program.  Figure  29  shows  the  result¬ 
ing  Us  tort  ions  compute  \. 

Approximate  the  extreme  near-wake  using  Equations  (2)  and  (3)  as 
shown  in  Figure  ~*8, 
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6.  Use  the  wake  from  Step  5  as  input  to  the  Prescribed  Wake  Program  and 
compute  a  second  estimate  of  tne  bound  circulation.  Compare  the  new 
circulation  distribution  with  the  approximation  used  in  Step  U  (see 
Figure  89).  If  differences  are  within  the  tolerances  indicated  in 
Appendix  II  to  be  acceptable  (as  is  the  case  in  Figure  89),  then  the 
wake  geometry  and  circulation  distribution  solutions  are  considered 
compatible;  if  not,  a  further  iteration  is  required  starting  with 
Step  3* 


EVALUATION  OF  ANALYTICAL  TECHNIQUES 


MODEL  BLADE  AIRFOIL  DATA 


The  airfoil  data  used  in  the  theoretical  calculation  were  based  on 
available  two-dimensional,  iow-Reynolds-number  data  adjusted,  through  a 
rynthesization  procedure,  to  provide  correlation  between  the  test  results 
for  the  untwisted  two-bladed  rotor  and  the  results  of  the  Blade-Element 
Momentum  Analysis.  Lerfornance  results  were  first  compared  on  Ct/ct  versus 
Cq/<r  plots  for  the  three  tip  Mach  numbers  tested  to  determine  the  general 
quality  of  correlation  between  test  and  this  particular  theory.  More 
detailed  comparisons  were  then  made  on  Cq>/a  versus  and  Cq/tr  versus 

8  plot.:  to  estimate  the  relative  magnitude  of  change  required  in  the 
airfoil  lift  curve  slope,  stall  angle,  and  drag  data  to  improve  correla¬ 
tion.  Approximately  five  iterations  were  required  at  each  tip  Mach 
number  condition  to  achieve  acceptable  correlation  in  both  the  thrust- 
torqae  results  and  the  thrust-torque-collective  pitch  results.  The  final 
synthesised  airfoil  data  generated  (Figure  90)  were  then  used  in  all  other 
theoretical  calculations.  In  comparing  these  data  with  the  original  two- 
dimensional  data,  only  small  changes  were  noted.  The  synthesized  data 
iiave  generally  higher  lift  curve  slopes,  higher  maximum  lift  coefficients, 
and  lower  profile  drag  coefficients.  The  differences  are  attributed  to 
differences  in  Reynolds  number  and  surface  conditions  for  the  blades 
tested. 


EVALUATION  <'Y  WAKE  GE'.fgTRV  PRCC.RAM 

•igures  U  through  95  compare  the  predicted  radial  and  ax;  al 
coordinates  of  the  “ip  vortex  with  those  measured  experimentally.  The 
experimental  curves  shown  were  obtained  from  the  faired  curves  given  in 
the  generalized  wage  charts  of  Figures  53  through  62.  The  basic  trends 
of  the  experimental  lata  have  already  been  discussed  in  the  section 
entitle*  V  .MtfS  -*v  EXPERIMENTAL  VTA  ICE  GEOMETRY  REAVITS ,  and  only  the 
general  ability  of  the  theory  to  predict  the  observed  results  is  considered 
here . 


Examination  of  the  results  in  Figures  91  through  95  leads  to  the 
following  general  observations  regarding  the  accuracy  of  the  theory: 

1.  '•'he  theory  predict s .values  of  k,  (the  average  slope  of  the  It  curve 

between  n  -  v...  -  “  .-r,  equivalently,  the  ncniiaens ional  axial 
ve 1  ~ci ty  which  are  significantly  lower  than  those  measured.  As  a 


result,  the  tip  vortex  from  one  blade  is  predicted  to  pass  closer 
to  the  following  blade  than  is  observed  experimentally.  This  is 
particularly  true  for  the  two-bladed  rotor  results,  where  the  pre¬ 
dicted  vortex-following  blade  distance  is  generally  about  one-half 
of  the  measured  value.  For  the  six-bladed  rotor  condition  shown  in 
Figure  93,  the  predicted  value  of  k^  is  actually  positive,  with  the 
result  that  the  vortex  passes  above  the  following  blade  instead  of 
below  the  blade  as  is  observed.  In  this  case,  the  predicted  and 
measured  vortex-to-following-blade  distances  are  about  the  same. 

2.  The  predicted  rate  of  contraction  of  the  tip  vortex  in  the  immediate 
vicinity  of  the  rotor  is  somewhat  less  than  that  measured.  As  a 
result,  the  predicted  radial  position  of  the  tip  vortex  at  the 
important  wake  azimuth  angle  of  360  deg/b  (i.e.,  in  tiie  vicinity  of 
the  following  blade )  is  approximately  0.02  R  farther  outboard. 

3.  The  analysis  predicted  accurately  the  experimental  values  of  ko,  the 
average  slope  of  the  z-j.  curve  for  ^>3^0  deg/b. 

4.  The  analysis  appears  to  be  generally  capable  of  predicting  the 
changes  observed  in  the  wake  geometry  when  various  rotor  or  flight 
condition  parameters  are  altered;  e.g.,  the  greater  k2  and  increased 
contraction  observed  with  increased  thrust  (Figures  91  and  93). 

The  results  obtained  from  the  Wake  Geometry  Program  for  the  full- 
scale  CH-h3A  rotor  were  also  compared  with  analytical  and  experimental 
results  presented  by  Clark  in  Reference  9.  The  comparison  is  shown  in 
Figure  96,  where  the  tip  vortex  axial  and  radial  coordinates  are  presented 
versus  wake  azimuth  angle.  Although  seme  difference  exists  in  the  rotor 
thru®*  levels,  the  comparisons  can  be  considered  quite  favorable,  at 
least  until  evidence  of  wake  instability  is  noted  in  the  results  of  this 
study  at  the  larger  values  of  >bv. 


In  summary,  the  wake  geometry  analysis  developed  herein  predicts 
many  of  the  qualitative  characteristics  observed  in  rotor  waxes.  The 
quantitative  accuracy  of  the  analysis  is  limited  by  one  major  shortcoming; 
namely,  the  inability  to  predict  accurately  the  average  axial  transport 
velocity  of  a  tip  vortex  element  between  the  time  when  it  is  generated 
and  the  time  when  the  element  passes  the  next  blade.  Although  the  exact 
reason  for  this  discrepancy  is  not  known  at  this  time,  the  need  for 
replacement  of  the  discrete  filament  wake  model  by  a  continuous  sheet 
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model  in  the  vicinity  of  the  blade  is  probably  implied.  Consideration 
should  also  be  given  to  reducing  the  computing  time  of  the  curre.it  Wake 
Geometry  Program  so  that  a  more  thorough  assessment  of  some  of  the  other 
assumptions  made  in  this  study  can  be  considered.  A  significant  reduction 
in  time  could  be  achieved  by  incorporating  the  symmetry  features  of  the 
hover  mode  of  operation. 


EVALUATION  OF  PERFORMANCE  METHODS 


The  theoretical  performance  methods  were  evaluated  by  comparing  the 
predicted  hover  performance  results  of  the  model  rotors  and  two  full- 
scale  rotors  with  test  results.  The  full-scale  rotors  selected  were  the 
six-bladed  Sikorsky  CH-53A  and  the  two-bladed  Bell  HU-LA.  The  discussion 
of  the  performance  results  follows. 

Model  Rotors 


Fterformance  Predicted  by  Uhcontracted  Wake  Methods 

The  performance  results  of  the  Momentum  Analysis,  Goldstein- Lock 
Analysis,  and  the  Prescribed  Classical  Wake  Analysis  are  presented  in 
Figures  97  through  102.  These  methods  do  not  account  for  wake  contraction 
and  will  hereafter  be  termed  uncontracted  wake  methods.  They  are  consid¬ 
ered  separately  frew  the  contracted  wake  methods  (Prescribed  Experimental 
Wake  Analysis  and  Prescribed  Theoretical  Wake  Analysis)  because  (l)  they 
represent  the  state-of-the-art  prior  to  the  consideration  of  contracted 
wake  effects,  ind  (2)  due  to  their  operational  simplicity  and  minimal 
computer  eost  requirements,  they  were  applied  to  a  greater  number  of  test 
conditions.  The  experimental  data  points  in  these  figures  represent  the 
results  of  the  fairings  of  the  collective  pitch  data  presented  in  Figures 
l5-’  through  18.  The  performance  results  for  varying  numbers  of  blades  are 
shown  in  Figures  97  and  for  the  -3-degree-twist ,  lB . 2-aspect-ratio 
rotor.  The  predicted  performance  ef  all  methods  is  in  good  agreement  with 
the  experimental  results  at  low  thrust  levels  but  becomes  increasingly 
optimistic  with  increasing  thrust  and  number  of  blades.  The  results  of 
the  three  undistortei  wake  methods  agree  closely  for  four,  six,  and  eight 
blades,  y 3r  two  blades,  some  separation  ef  the  predicted  results  at  the 
high  thrust  levels  is  evident.  The  aforementioned  results  were  found  to 
be  generally  true  for  all  the  rotors  tested,  ana  the  trends  with  thrust 
level  and  number  of  blades  are  the  same  as  have  been  cited  for  full-scale 
r 'tors  dr.  Reference  1.  Correlation  results  for  rotors  of  C  and  -16 
degrees  of  twist  are  presented  in  Figures  >'•>  and  IOC,  respectively,  for 


two  and  six  blades.  It  should  be  noted  that,  as  discussed  previously,  the 
airfoil  data  used  was  synthesized  by  matching  the  results  of  the  Momentum 
Analysis  to  the  test  results  for  the  two-bladed,  zero  twist  rotor.  For 
six  blades,  the  predicted  performance  is  increasingly  optimistic  with 
decreasing  blade  twist  (e.g,,  compare  results  from  Figures  99  and  100  at 
Ct/ct  -  0.075)*  In  Figure  101,  the  results  for  the  low-aspect-ratio 
blades  (AR  =  13*6)  are  presented.  Slightly  improved  correlation  between 
theoretical  and  experimental  results  compared  to  that  of  the  high-aspect- 
ratio  blades  (AR  =  18.2)  of  Figure  97  is  indicated.  Since  fo’-  the  same 
number  of  blades  and  rotor  radius,  rotor  solidity  increases  with  decreas¬ 
ing  aspect  ratio,  this  indicates  that  the  correlation  improves  with 
increased  rotor  solidity  when  such  increases  result  from  varying  aspect 
ratio.  This  is  opposite  to  the  correlation  trend  established  when 
solidity  was  increased  by  varying  the  number  of  blades.  The  dependence 
of  the  correlation  on  the  manner  in  which  solidity  is  varied  (i.e.,  number 
of  blades  versus  aspect  ratio)  also  agrees  with  the  conclusions  of 
Reference  1  based  on  full-scale  rotors.  Finally,  the  small  improvement 
in  correlation  of  the  predicted  and  measured  model  performance  results  at 
a  tip  speed  of  52c-  fps  (Figure  102),  noted  for  the  two-bladed  rotor,  also 
parallels  the  full-scale  results  of  Reference  1,  in  which  tho  correlation 
is  shown  to  improve  with  decreasing  tip  Mach  number. 

It  is  thus  concluded  that  the  discrepancies  between  the  experimental 
performance  and  the  performance  predicted  by  the  uneontracted  wake  methods 
generally  exhibit  the  same  trends  for  the  model  rotor  as  they  do  for  full- 
scale  rotors:  the  deterioration  in  performance  prediction  with  in  Teasing 
number  of  blades,  blade  loading,  and  tip  Mach  number.  In  addition,  a 
tendency  for  improved  agreement  with  increased  blade  twist  was  noted  with 
the  models.  These  discrepancies  are  believed  to  be  primarily  the  result 
of  the  assumptions  made  in  these  methods  regarding  the  rotor  wake  geoeetry 
and  the  impact  which  such  assumptions  have  on  the  proximity  of  the 
operating  conditions  of  the  tip  sections  to  stall. 

Integrated  Performance  Predicted  by  the  Contracted  Wake  Analyses 

The  performance  results  from  Uje  contracted  wake  methods  (i.e.,  the 
Prescribed  Experimental  Wake  Analysis  and  the  Prescribed  Theoretical  Wake 
Analysis)  are  compared  with  the  test  results  in  Figures  103  through  106. 
Also  included  for  comparison  are  the  results  from  one  of  the  uncontracted 
wake  analyses  ((Ooldstein-Loek  Analysis)  presented  previously.  The 
Goldstein-Lock  Analysis  was  selected  because  it  is  a  widely  used  hover 
performance  analysis.  The  results  for  all  conditions  analyzed  using  the 
contracted  wake  methods  are  included  in  these  figures.  Application  of  the 
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Fteccribed  Theoretical  Wake  Analysis  to  a  six-bladed  model  rotor  was 
limited  to  one  condition  because  of  the  computing  time  requirements  of 
the  program  at  this  time.  For  consistency,  the  generalized  experimental 
wake  results  of  Figures  53  through  65  were  used  to  obtain  all  Presciibed 
Experimental  Wake  method  results  shown  in  this  set  of  figures.  The 
sensitivity  of  predicted  performance  to  variations  from  the  generalized 
wake  coordinates  was  significant  and  will  be  discussed  later  in  this 
section.  The  performance  predicted  by  the  prescribed  contracted  wake 
methods  generally  correlates  well  with  the  experimental  model  rotor 
results,  particularly  for  the  six-bladed  rotors.  Reasonable  correlation 
with  collective  pitcn  also  exists  except  for  the  low-aspect-ratio  blades 
(Figure  105 ),  for  which  all  of  the  theories  predict  higher  thrust  and 
torque  than  measured  for  a  given  pitch.  Except  in  one  instance,  the 
thrust  based  on  the  theoretical  wake  is  consistently  lower  than  that  pre¬ 
dicted  using  the  experimental  wake.  This  is  due  to  the  previously 
mentioned  prediction  of  the  tip  vortex  as  being  nearer  the  blade  than  was 
observed  from  experiment.  The  prescribed  contracted  wake  methods  appear 
to  have  one  distinct  advantage  over  the  uncontracted  wake  methods :  the 
ability  to  predict  more  accurately  the  performance  trends  with  increased 
number  of  blades  and  blade  loading  (C.T./cr).  As  shown  previously  in 
Figures  97  through  102,  the  Momentum,  Goldstein-Lock,  and  Prescribed 
Classical  Wake  methods  all  predict  optimistic  performance  for  high  thrust 
'  vels  and  high  number  of  blades.  For  the  six-bladed  rotor  operating  at 
the  high  thrust  condition,  the  results  of  the  Prescribed  Experimental 
Wake  Analysis  based  on  the  generalized  wake  indicate  reduced  performance 
relative  to  the  Goldstein-Lock  results  and  improved  correlation  with  the 
measured  data. 

Section  Characteristics 


Some,  insight  into  the  reasons  for  the  above-mentioned  changes  in 
predicted  performance  may  be  gained  from  Figure  107,  where  the  theoretical 
spanvise  distributions  of  axial  induced  velocity,  section  angle  of  attack, 
thrust,  and  torque  are  compared.  The  predicted  distributions  correspond¬ 
ing  to  the  six-bladed  rotor  of  Figures  97  and  103  operating  at  a  10-degree 
collective  pitch  setting  were  selected  for  this  comparison.  The  influence 
of  the  experimental  tip  vortex  geometry  on  the  induced  velocity  distribu¬ 
tion  and  the  resulting  effects  on  the  angle  of  attack,  thrust,  and  torque 
distributions  are  apparent  in  Figure  107-  For  lew  thrust  levels,  number 
of  blades,  and  tip  Mach  numbers ,  the  distributions  of  the  uncoiutractfe4 
wake  methods  are  generally  compensative  between  the  inboard  and  tip 
regions  of  the  blades  as  far  as  integrated  performance  is  concerned. 
However,  at  a  condition  such  as  i.h*t  of  this  figure,  the  contracted  wake 
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geooetry  considerations  result  in  a  decrease  in  integrated  performance 
relative  to  that  predicted  by  the  uncontracted  wake  methods.  In  the  real 
wake,  the  tip  vortices  are  positioned  much  closer  to  the  tip  path  plane 
than  is  assumed  in  the  Goldstein-Lock  and  Prescribed  Classical  Wake 
analyses.  As  the  number  of  blades  increases,  due  to  the  reduced  separa¬ 
tion  between  blades,  each  blade  is  closer  to  the  tip  vortex  generated  by 
the  blade  ahead,  and  increased  aerodynamic  interference  results.  The 
circulation  strength  of  the  tip  vortex  also  increases  with  the  increased 
thrust  level  associated  with  high  manbers  of  blades.  This,  coupled  with 
the  close  proximity  of  the  preceding  blade's  oip  vortex  to  the  blade, 
results  in  increased  local  angle- of-attack  values  in  the  blade  tip  region 
due  to  the  upflow  generated  by  the  contracted  vortex.  It  should  be  noted 
that  the  predicted  performance  deterioration  is  attributable  to  deterio¬ 
ration  in  integrated  thrust  rather  than  torque.  In  fact,  the  predicted 
torque  based  on  the  experimental  wake  for  the  subject  test  condition  is 
slightly  less  than  the  torque  predicted  by  the  undistorted  wake  analyses. 
As  shown  in  Figure  107c,  this  is  caused  by  the  reduction  in  induced 
torque. 

ferformance  Sensitivity  to  Tip  Vortex  Geometry 

As  mentioned  previously,  the  performance  results  were  found  to  be 
sensitive  to  certain  wake  parameters.  The  sensitivity  results  are 
summarized  in  the  following  listing,  in  which  the  wake  parameters  and  the 
their  primary  influence  on  the  wake  geooetry  ere  presented  in  decrensing 
order  of  importance : 

k^  and  k2  (the  axial  position  of  the  tip  vortex) 

A  and  A  (the  radial  position  of  the  tip  vortex) 

Klf  =  i  110(1  K2y  -  i  (the  1x1111  position  of  the  outer  region  of 

the  vortex  sheet) 

Ki_  „  and  K0_  (the  axial  position  of  the  inner  region  of 

=  O  =  0 

the  vortex  sheet) 

The  axial  position  of  the  tip  vertex,  defined  by  k^  ar*d  k2,  was  found  to 
have  the  dominant  influence  on  the  performance  results.  To  illustrate  the 
sensitivity  to  the  tip  vortex  axial  position,  the  effect  of  variations  in 
on  hover  performance  is  shewn  in  Figure  1^.  k^  was  varied  so  as  to 
axially  displace  the  tip  vortex  from  the  generalized  wake  position  by 
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-Vfs  of  the  rotor  radius  (-l'jt  R)  at  and  beyond  the  wake  azimuth  position 
equal  to  the  blade  spacing4.  This  effectively  moved  the  tip  vortex  1%  R 
closer  to  or  farther  from  the  following  blade.  The  magnitude  of  the 
changes  approximately  represents  the  accuracy  of  the  tip  vortex  general¬ 
isation.  The  performance  results  show  increasing  sensitivity  with  number 
of  blades  (two  versus  six  in  Figure  103).  To  analyze  these  results  in 
greater  detail,  distributions  of  the  blade  section  characteristics  were 
compared.  For  example,  the  blade  airload  distribution  and  tip  vortex 
positions  are  shown  in  Figure  109  for  the  lettered  conditions  indicated 
on  Figure  10S.  (Note  that  the  nominal  conditions  A-D  are  the  ones  for 
which  tip  vortex  streamlines  were  presented  in  Figure  ?0.)  The  larger 
change  in  loading  distribution  for  the  six-bladed  rotor  conditions  (C  and 
C"  ;  D  and  D"  )  is  due  to  the  increased  proximity  of  the  tip  vortex  to  the 
blade.  Also  evident  in  this  figure  is  (l)  the  typical  change  in  the 
character  of  the  loading  in  the  tip  region  with  varying  number  of  blades 
and  thrust  level,  and  (2)  the  radially  outward  movement  of  the  peak 
loading  with  increased  number  of  blades  which  corresponds  to  the  movement 
in  the  radial  position  of  the  tip  vortex. 

Predicted  performance  was  fairly  insensitive  to  changes  in  A  and  A 
providing  they  were  changed  in  combination  rather  than  independently  so 
as  to  maintain  the  radial  coordinates  of  the  near-wake  portion  of  the  tip 
vortex  within  the  experimental  accuracy  --  particularly  beneath  the 
following  blade.  As  mentioned  previously,  A  =  0.78  was  selected  in  the 
wake  generalization  to  curve-fit  the  near-wake  radial  coordinates  when 
used  in  combination  with  the  specific  A  values  given  in  Figure  62  as  a 
function  of  Cx*  For  a  given  test  condition,  there  are  other  combinations 
of  A  and  A  that  would  fit  the  near-wake  data  within  the  experimental 
accuracy.  However,  the  performance  was  relatively  insensitive  to  such 
changes  even  though  the  far-wake  contraction  varied  considerably  (0.78  to 
0.72).  least  senstivity  resulted  from  variations  of  the  inboard  wake 
(vortex  sheet).  For  example,  a  10%  change  in  the  axial  coordinates  of 
the  vortex  sheet  for  a  typieal  condition  resulted  in  less  than  a  1%  change 
in  thrust  and  torqus  for  a  two-bladed  rotor  and  less  than  a  2%  change  for 


•Results  for  a  1%  R  displacement  of  the  vip  vortex  of  the  sir-bladed  rotor 
toward  the  rotor  disc  are  not  included  in  Figure  100.  This  displacement 
resulted  in  the  vortex  passing  within  0.^4  R  of  the  following  blade. 
Theoretical  applications  in."’l ring  passage  of  the  vortex  within  this 
distance  cf  the  blnde  are  questionable  due  to  the  limitations  of  lifting 
line  theory  and  limited  information  or.  vortex  core  site. 


a  six-bladed  rotor.  In  summary,  rotor  performance  becomes  more  sensitive 
to  wake  geometry  variations  as  number  of  blades  and/or  thrust  level  is 
increased.  Performance  is  particularly  sensitive  to  the  axial  position 
of  the  tip  vortices  due  to  their  close  proximity  to  the  blades.  This 
sensitivity  of  blade  loading  to  small  changes  in  tip  vortex  position 
predicted  by  the  Prescribed  Wake  Analysis  confirms  certain  phenomena 
observed  on  full-scale  rotors.  This  will  be  discussed  in  the  following 
section. 

Full-Scale  Rotors 


The  theoretical  methods  were  applied  to  compute  the  hover  performance 
of  the  six-bladed,  Sikorsky  CH-53A  rotor.  Two-dimensional  airfoil  data, 
supplied  by  Sikorsky,  were  used  in  all  calculations.  In  Figure  110,  the 
theoretical  performance  results  are  compared  with  test  results  from  the 
Sikorsky  whirl  stand  shown  in  Figure  7.  The  absolute  accuracy  of  thrust 
at  a  given  power  level  measured  on  this  facility  is  estimated  as  -2^. 

The  test  data  presented  in  Figure  110  have  been  corrected  for  ground 
effect  and  whirl  stand  interference.  As  shown  in  this  figure,  the  perfor¬ 
mance  predicted  by  the  Momentum,  Gollstein-Lock,  and  Prescribed  Classical 
Wake  analyses  (uncontractei  wake  methods)  js  increasingly  optimistic  with 
increasing  thrust  level,  which  is  consistent  with  the  previously  recog¬ 
nized  trend.  On  the  other  hand,  the  results  of  the  Prescribed 
Experimental  Wake  and  Prescribed  Theoretical  Wake  analyses  (contracted 
wake  methods)  are  slightly  pessimistic  relative  to  the  experimental  data. 
The  generalized  wake  coordinates  from  the  model  rotor  test  were  used  in 
the  Prescribed  Experimental  Wake  Analysis  to  obtain  the  full-scale  CH-53A 
performance  results.  The  performance  corresponding  to  small  deviations 
of  the  tip  vortex  geometry  from  the  generalized  coordinates,  within  the 
experimental  wake  accurccy,  is  presented  in  Figure  111.  It  is  shown  that 
the  performance  veriations  associated  with  small  changes  in  wake  geometry 
are  significant.  In  fact,  as  indicated,  it  was  possible  to  obtain 
correlation  with  the  test  data  by  displacing  the  tip  vortex  only  R 
farther  away  from  the  rotor  disc. 

The  section  angle  of  attack  distributions  predicted  by  the  various 
theoretical  methods  for  the  12-degree  collective  pitch  condition  of 
Figure  110  are  presented  in  Figure  112.  The  prescribed  experimental  wake 
results  again  exhibit  a  large  angle-of-attack  increase  near  the  tip 
caused  by  the  strong  influence  of  the  contracted  tip  vortex.  The  angle- 
of-attack  distribution  in  the  tip  region  of  the  blade  has  been  substan¬ 
tiated  on  the  Sikorsky  whirl  stand  in  a  recent  test  in  which  pressure 
measurements  were  recorded  at  three  stations  over  tie  outer  of  the 
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blade.  It  was  found  that  with  minor  adjustments  to  the  generalized  wake 
coordinates,  the  angle-of-attack  distribution  and  total  rotor  performance 
could  be  predicted  within  the  experimental  limits  of  the  test  data.  This 
emphasizes  the  inadequacy  of  *he  uncontracted  wake  methods  for  predicting 
the  distributions  of  the  performance  characteristics.  Also,  although  this 
is  encouraging  insofar  as  the  accuracy  of  the  Prescribed  Experimental 
Wake  Analysis  is  concerned,  additional  pressure  measurement  data  suffi¬ 
cient  in  extent  to  obtain  the  complete  angle-of-attack  distribution  along 
the  blade  for  various  rotors  and  operating  conditions  would  obviously  be 
desirable  to  permit  further  evaluation  of  this  method  and  refinement  of 
the  generalized  wake  model.  Examination  of  the  theoretical  radial  loading 
distributions  associated  with  the  operating  conditions  of  Figure  112 
revealed  that  the  outer  1056  of  the  blade  generates  approximately  one-third 
of  the  blade  lift  according  to  this  theory,  compared  to  only  one-quarter 
of  the  lift  according  to  the  Goldstein-Lock  theory.  This  comparison 
implies  an  extreme  sensitivity  of  rotor  thrust  to  the  aerodynamic  charac¬ 
teristics  over  the  blade  tip  region  and  the  need  for  an  accurate  simula¬ 
tion  of  the  airfoil  characteristics  and  velocities  induced  by  the  waku  in 
this  area.  The  assumption  of  two-dimensional  flow  over  the  blade  span 
has  been  shown  in  Reference  2k  to  be  questionable  in  the  blade  tip  region. 
The  magnitude  cf  three-dimensional  corrections  and  their  influence  on  the 
results  predicted  by  the  prescribed  wake  analyses  remain  to  be  determined. 

The  theoretical  methods  were  also  applied  to  compute  the  hover 
performance  of  the  two-bladed  Bell  HU-LA  rotor.  Again,  the  Blade  Element- 
Momentum,  Goldstein-Lock,  and  Prescribed  Wake  methods  were  used.  Two- 
dimensional  airfoil  data,  synthesized  by  the  Bell  Helicopter  Company,  were 
used  for  the  NACA  0015  airfoil  section  in  all  calculations.  In  Figure 
113,  the  theoretical  performance  results  are  compared  with  test  data  from 
a  tethered  aircraft  test  reported  in  Reference  25  (whirl-stand  data  for 
isolated  rotor  performance  was  not  available  for  two-bladed  rotors).  For 
the  test  data  shown,  the  engine  shaft  horsepower  was  reduced  by  15%  to 
account  for  tail  rotor,  transmission,  and  accessory  power  losses.  The 
wake  lata  used  in  the  prescribed  experimental  wake  calculations  were 
obtained  by  interpolating  the  generalized  rotor  wake  data  to  the  appro¬ 
priate  twist  and  thrust  levels.  The  theoretical  results,  including 
contracted  wake  results,  are  all  optimistic  relative  to  the  experimental 
results.  It  is  somewhat  uncertain,  for  this  rotor,  as  to  whether  the 
performance  discrepancies  are  attributable  to  limitations  of  the-  theoret¬ 
ical  methods,  to  inaccuracies  in  the  estimation  of  the  power  losses, 
and  ’or  to  the  synthesis  of  the  airfoil  data.  It  is  thus  suggested  that 
the  results  presented  in  Figure  113  be  considered  preliminary  until  other 
two-feladed  rotors  are  analysed. 


The  sensitivity  of  blade  loading  to  small  changes  in  the  tip  vortex 
position  tends  to  confirm  vortex  interferenc'  as  the  source  of  certain 
phenomena  observed  on  full-scale  rotors.  As  explained  in  Reference  1, 
the  presence  of  a  small  amount  of  ambient  wind  can  produce  changes  in  the 
tip  vortex  position  which,  at  certain  conditions,  are  reflected  in  the 
rotor  characteristics  as  vibratory  flapping  and  tip  stall.  Since  the 
publication  of  that  reference,  blade  tracking  problems  at  high  thrust 
levels  (greater  than  normal  operating  levels)  ha/e  been  related  to 
vortex  interference  effects.  The  Prescribed  Wake  Program  appears  Lo  be 
of  potential  value  for  analyzing  these  phenomena.  However,  for  increasing 
number  of  blades,  the  problem  associated  with  the  sensitivity  of  the 
predicted  results  to  minor  wake  geometry  variations  is  recogni'  ?d.  With 
further  experience  with  the  method,  and  through  small,  systematically 
developed  adjustments  to  the  generalized  wake  model  presented  herein 
and/or  further  refinements  in  the  Wake  Geometry  Program,  it  may  become 
possible  to  consistently  predict  hover  performance  to  a  high  degree  of 
accuracy  with  the  method  in  its  present  form.  However,  until  this  is 
accomplished,  further  investigation  of  three-dimensional  tip  effects  and 
unsteady  wake  effects  associated  with  both  ambient  wind  variations  and 
the  apparent  instability  of  ti.c  far  wake  appears  warranted. 


RESULTS  AND  CONCLUSIONS 


1.  Model  rotors  of  the  scale  tested  can  be  successfully  used  to 
obtain  performance  and  wake  geometry  data  which  are  both  systematic  and, 
more  importantly,  indicative  of  those  characteristics  observed  for  full- 
scale  rotors.  The  use  of  an  indoor,  small-scale  facility  greatly  reduces 
the  cost  of  acquiring  such  information. 

2.  The  model  rotor  wake  geometry  for  the  region  directly  beneath 
the  rotor  (i.e.,  within  approximately  one-fourth  of  the  rctor  radius)  can 
be  represented  by  simple  generalized  equations  which  facilitate  the  rapid 
estimation  of  rotor  wake  coordinates  for  a  wide  range  of  rotor  designs  and 
operating  conditions.  The  generalization  was  based  on  the  following 
observations  : 

a.  The  rate  of  descent  of  an  element  of  the  tip  vortex  from  a 
blade  is  substantially  constant  prior  to  its  passage  beneath 

the  following  blade.  This  axial  velocity  increases  approximately 
linearly  with  increasing  blade  loading  and  decreases  slightly 
with  increasing  blade  twist.  Within  the  experimental  accuracy 
of  the  data,  the  axial  velocity  is  insensitive  to  variations  in 
number  of  blades. 

b.  The  axial  velocity  of  a  tip  vortex  element  after  the  passage  of 
the  following  blade  is  increased  to  a  new,  relatively  constant 
value.  This  velocity  increases  with  increasing  disc  loading  in 
a  manner  proportional  to  the  momentum  inflow  velocity.  The 
constant  of  proportionality  is  approximately  l.U  for  rotors  with 
untwisted  blades  and  decreases  with  increasing  twist. 

c.  Hte  radial  position  of  a  tip  vortex  element  decreases  in  a 
decaying  exponential  manner.  The  rat#  of  decrease  (or  contrac¬ 
tion)  appears  to  be  primarily  determined  by  rotor  disc  loading 
and  is  greater  with  increasing  disc  loading. 

i.  The  inboard  vortex  sheet  from  each  blade  is  fairly  linear  with 
radial  position  at  a  specific  wake  azimuth  location.  The  axial 
velocity  of  the  sheet  increases  with  increasing  disc  loading  and, 
to  a  lesser  extent,  with  blade  twist.  The  rate  of  descent  of 
the  outer  portion  of  the  vortex  sheet  is  approximately  twice 
that  observed  for  the  tip  vortex  following  *hc  passage  of  the 
next  blade. 


e.  Wake  geometry  is  insensitive  to  independent  variations  in  aspect 
ratio  and  tip  speed.  The  influence  cf  number  of  blades  is 
limited  to  the  establishment  of  the  wake  azimuth  angle  at  which 
the  wake  axial  descent  velocities  are  observed  to  increase 
significantly. 

3.  The  stability  of  the  experimental  wake  decreases  with  increasing 
distance  from  the  rotor,  indicating  that  the  classical  concept  of  a 
smoothly  contracting  wake  may  be  incorrect. 

4.  Attempts  to  develop  an  analytical  method  for  predicting  the 
contracted  tip  vortex  geometry  based  upon  the  interaction  of  discrete 
vortex  filaments  were  partially  successful.  The  method  predicts  the 
general  features  of  the  tip  vortex  and  the  qualitative  variation  of  these 
features  with  changes  in  rotor  or  flight  condition  parameters,  however, 
the  position  of  the  tip  vortex  shed  by  one  blade  relative  to  the  following 
blade  was  not  accurately  predicted  by  the  analysis.  Compared  to  the 
experimental  observations,  the  predicted  vortex  was  generally  closer  to 
the  blade  and  farther  outboard.  In  addition,  the  contracting  tip  vortex 
was  predicted  to  become  unstable  at  moderate  dis.-nces  below  the  rotor 
plane ,  thus  complicating  the  problem  of  predicting  the  wake  geometry 
beyond  the  near-wake  region.  Available  experimental  evidence  appears  to 
substantiate  this  prediction. 

5.  Rotor  hover  performance  methods  which  assume  that  the  blades  are 
lifting  lines  and  that  the  rotor  wake  is  uncontracted  give  generally 
comparable  performance  predictions.  As  was  noted  for  full-scale  rotors, 
the  accuracy  of  these  methods  in  predicting  model  rotor  performance 
decreased  as  the  number  of  blades,  blade  loading,  and  tip  Mach  number 
increased.  In  addition,  the  model  rotor  correlations  indicated  a 
decrease  in  accuracy  when  the  twist  of  the  rotor  blades  under  consideration 
was  reduced.  The  partial  success  achieved  with  these  methods  in  the  past 
results  from  compensating  errors  in  the  blade  loading  distributions. 

6.  The  elimination  of  the  uncentracted  wake  assumption  through  the 
use  of  generalized  experimental  wake  data  significantly  improves  the 
accuracy  of  performance  characteristics  predicted  for  model  rotor  condi¬ 
tions  for  which  the  uncontracted  wake  methods  exhibit  major  shortcomings 
(i.e.,  for  conditions  with  increasing  number  of  blades  and  thrust  levels). 
Application  of  the  experimental  wake  method  to  full-scale  rotors,  using 
model  rotor  wake  data,  yielded  similar  improvements  for  the  limited 
number  of  conditions  examined.  In  addition,  the  method  predicts  blade 
spanwise  angle-of-attack  distributions  which  are  greatly  different  and 
more  realistic  than  those  predicted  using  uncontracted  wake  analyses. 
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7.  The  use  of  predicted  contracted  wake  geometries  in  computing 
rotc-r  performance  led  generally  to  reasonable  but  more  conservative 
(lower  thrust  for  same  power)  estimates  relative  to  those  obtained  using 
the  generalized  experimental  wake  results.  This  was  mainly  due  to  the 
prediction  of  a  smaller  axial  distance  between  the  tip  vortex  and  the 
following  blade  than  was  observed  experimentally. 

3.  The  performance  predicted  on  the  basis  of  lifting  line  theory 
using  a  realistic  contracted  wake  is  very  sensitive  to  small  changes  in 
wake  geometry,  particularly  for  rotors  for  which  the  tip  vortex  is 
positioned  very  close  to  the  following  blade.  While  the  generalized  wake 
geometry  results  presented  herein  are  believed  to  be  based  on  reasonable 
fairings  of  the  data  available,  in  view  of  this  sensitivity,  small 
systematic  refinements  to  the  fairings  may  be  required  if  consistently 
accurate  performance  predictions  are  to  be  achieved. 


RLCCMMENDATIONS 


1.  Because  of  the  more  accurate  definition  of  the  radial  distribu¬ 
tions  of  blade  section  angles  of  attack  which  they  provide,  contracted 
wake  analyses  should  be  used  to  predict  hover  performance.  They  should  be 
particularly  useful  in  evaluating  and  developing  improved  blade  tip 
designs. 

2.  In  lieu  of  an  accurate  method  for  predic  .  ..g  the  geometry  of  a 
contracted  wake,  the  Prescribed  Wake  Analysis,  together  with  the 
generalised  experimental  wake  geometry  charts,  should  be  employed. 

3.  Efforts  to  obtain  accurate  wake  geometry  data  from  full-scale 
rotors  should  continue  to  substantiate  further  the  applicability  of 
generalized  model  rotor  wake  geometry  data  to  full-scale  rotors. 

4.  Measurements  of  blade  pressure  distributions  should  be  obtained 
to  provide  experimental  data  for  detailed  comparisons  with  predicted 
distributions  based  on  the  generalized  wake  information. 

5.  An  experimental  investigation,  employing  model  rotor  flow 
visualization  techniques,  should  be  conducted  to  obtain  systematic  rotor 
performance  and  wake  geometry  data  for  blade  designs  with  taper,  nonlinear 
twist,  and  premising  tip  shapes. 

6.  Investigations  should  be  undertaken  to  examine  in  detail  the 
stability  characteristics  of  the  wake  of  a  hovering  rotor. 

7.  Further  work  should  be  undertaken  to  reduce  the  computing  time 
required  by  the  analysis  developed  herein  for  predicting  the  geometry 

of  the  rotor  wake.  This  will  greatly  facilitate  further  studies  designed 
to  improve  the  quantitative  accuracy  of  \he  analysis. 
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Figure  3-  Schematic  Cress  Section  of  Kotor  Te§+.  Rig 
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Figure  4.  Rotor  Hub. 
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Schematic  of  Model  Rotor  Blade  Construction 


Figure  6.  Flow  Visualization  Grid 
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Figure  12.  Tip  Vort e*  Coordinates  Measured  From  Photos  of  Figure  11 
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Figure  15 .  Experimental  Model  Kotor  Performance  —  9 
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Figure  15-  Continued 
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Figure  15-  Concluded. 
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(a)  HR  -  TOO  FF5 

Figure  16.  Experimental  Model  Rotor  Rgrforaanee  — 
=  -8°,  AR  =  18.2. 
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Figure  IT-  Concluded. 
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(a)  flR  =  700  FF5 

Figure  18.  Experimental  Model  Rotor  Performance  — 
§  x  =  -8*  ,  AR  =  13.6. 
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TORQUE  COEFFICIENT  SOLIDITY,  C» 

(b)  11 R  =  525  TVS 


Figure  i8.  Concluded. 
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Figure  19.  Typical  Effect  of  Collective  Pitch  or;  Model  Rotor  Perfornance . 
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£ur*  ?0.  Typical  Experimental  hodel  Kotor  Perform- cc  Expressed 
ir.  Teres  of  Thrust  and  Torque  Coefficients. 
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Figure  23.  Effect  of  Tip  Speed  on  Experimental  Model  Rotor  Performance 

^  *  0*  ,  AR  -  18.2. 
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Figure  21*.  Effect  of  Tip  Speed  011  Experimental  Model  Rotor  Performance  - 

S  s  -l6‘,  AR  =  18.2. 
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PLASTIC  TIP  BLADES 


perimental  Model  Rotor  Performance  for  Various  Fetor  Heights  Above  the  Ground. 


Effect  of  Rotor  Height  Above  the  Ground  on  Rotor  Thrust 
Augmentation . 
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Figure  27.  Sequence  of  ftictographu  Showing  Time  History  of  V'ake  --  b  «  1* . 
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Figure  30.  Sequence  of  Pantographs  Showing  Wake  Instability. 
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(b)  *  =  90* to  120* 
Figure  30.  Concluded. 
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Figure  32 .  Typical  Wake  Flow  Time  Exposure  Photograph, 
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Figure  3^.  Photographs  Showing  Effect  of  Number  of  Blades 
and  Thrust  Level  on  Wake  Geometry. 
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Figure  35.  Photographs  Showing  Effect  of  Number  of  Blades  on 
Wake  Geometry  for  Rotor  Solidity  of  O.lUo. 
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Photograph  Showing  Effect  of  Ground  and  Whirl  Stan< 
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QR  700  FP;  C  T  ~  0.0034 
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Figure  1+5 .  Effect  of  Number  of  Blades  on  Tip  Vortex  Coordinates  for  Constant  Disc  Loading 
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Iffect  of  Twist  on  Tip  Vortex  Coordinates 


Figure  1*9.  Effect  of  Tip  Speed  on  Tip  Vortex  Coordinates 
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Figure  U9.  Concluded. 
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QR  =  700  F  PS  C,/cr~0.10  CT~0.0035 
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Figure  51.  Effect  of  Whirl  Stand  Conditions  on  Tip  Vortex  Coordinates 
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INBOARD  WAKE  PARAMETER 
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Figure  6U.  Experimental  ^  Inboard  Wake  Parameter. 
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"igure  *7.  Comparison  of  'Generalised  Wake  Results  With  Those 
of  References  1,  10,  and  11  --  Axial  Coordinates. 
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Figure  68.  Ccwparison  of  Generalized  Wake  Results  With  Those 
of  References  1,  10,  ar.d  11  --  Radial  Coordinates. 
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Figure  70.  Comparison  of  Generalized  Wake  Boundaries 
With  Results  of  Reference  l4. 
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Figure  71.  LT<\RL  Prescribed  Wake  Hover  Performance  Program. 


152 
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Figure  73 •  Iteration  Procedures  for  Computing  Wake  Geometries. 
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Figure  74.  Segmented  Discrete  Vortex  Representation  of  the  Wake 
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Figure  75.  Typical  Truncated  Wake  Results. 
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Figure  77.  Coordinate  Systems  and  Nomenclature  for 
Theoretical  Wake  Geometry  Computations. 
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Figure  78.  Intersection  of  Wake  With  r-z  Plane  Shoving 
Typical  Inboard  Vortex  Sheet  Representation, 
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30  rotation  of  rotor 


TIP  VORTEX  RADIAL  COORDINATE,? 

Figure  79.  Typical  Predicted  Rotor  WaJce  Boundaries  At  Various  Time  Steps 


Figure  80.  Model  Rotor  Photographs  Sheering  l^rpical 
Local  ^oil-up  of  Tip  Vortex. 
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NOTE 

ROLL-UP 


Figure  8l.  Photograph  of  Wake  for  a  Full-Scale  Rotor 
Showing  Local  Roll-Up  of  the  Tip  Vortex. 
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Figure  82.  Photograph  of  Wake  of  Model  Rotor  Taken  in  a  Water  Tunnel 
Shewing  Local  Roll-Up  of  the  Tip  Vortex. 


165 


166 


STABLE 


0 


50 


10  20  30  40 

TIME  STEP  (1  TIME  STEP  =  30°  ROTATION  OF  ROTOR) 

Figure  84.  Typical  Time  Histories  of  Computed  Tip  Vortex  Coordinates 
Showing  Convergence  of  Extreme  Near-Wake. 
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TIP  VORTEX  COORDINATES,  r  AND 


EXPERIMENT 


WAKE  AZIMUTH  COORDINATE,  <//w  r'EG 


Figure  85.  Computed  Radial  and  Axial  Coordinates  of 
Tip  Vortex  for  Two  Time  Steps. 
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Figure  87.  Approximation  to  Blade  Circulation 


Distribution  in  Sample  Computation. 
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Figure  88.  Approximation  to  Extreme  Near-Wake  in  Sample  Computation. 
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Figure  91.  Predicted  Effect  of  Rotor  Thrust  on  Wake  Geometry. 
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Figure  9?.  Predicted  Effect  of  31ade  Twist  on  Wake  Geometry. 
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Figure  93.  Predicted  Effect  of  Number  of  Blades  and  Thrust 
Level  on  Wake  Geometry. 
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Figure  95.  Predicted  Effect  of  Tip  Speed  on  Wake  Geometry. 
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Figure  96.  Comparison  of  Predicted  Wake  Geometry- 
Results  With  Results  of  Reference  9. 
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figure  97.  Comparison  of  Results  of  Unccntracted  Wake  Analyses  With 
Experimental  Performance  Results  for  Two-  and  Six-Bladed 
Rotors  --  =  -8°,  AR  =  l8.2,  HR  =  700  fps. 
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Figure  90*  Comparison  of  Results  of  Uncontracted  Wake  Analyses  With 

Experimental  Performance  Results  for  Four-  and  Eight-Bladed 
Model  Motors  --  I?  =  -8°,  AR  =  .?,OR  =  700  fps. 
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Figure  100.  Coaparison  of  Results  of  Uhcortraeted  Wake  Analyses  With 


Experinental  Performance  Results  for  Two-  and  Six-Bladed 
Model  Rotors  —  0^  *  -16°,  AR  =  18.2,  fiR  =  TOO  fps . 
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Figure  101.  Comparison  of  Results  of  Uncontracted  Wake  Analyses  With 
Experimental  Performance  Results  for  Two-  and  Six-Bladed 
Model  Rotors  —  =  -8°,  AR  =  13.6, fiR  =  700  fps. 
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Experimental  Performance  Results  for  Two-  and  Six-Bladed 
Model  Rotors  —  Q  =  -8°,  AR  =  i8.2,  HR  =  525  fps . 
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Figure  103.  Comparison  of  Results  of  Contracted  Wake  Analyses  >Iith 

Goldstein-Lock  Results  and  Experimental  Performance  Results 
for  Model  Rotors  —  0^  =  -8°,  AR  =  18.2,  OR  =  700  fps . 
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THRUST  COEFFICIENT/SOLIDITY,  CT/ 


Figure  104.  Comparison  of  Results  of  Contracted  Wake  Analyses  With 

Goldstein-Lock  Results  and  Experimental  Performance  Results 
for  Model  Rotors  --  6^  =0°,  AR  =  18.2,  HR  =  700  fps. 
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Figure  105.  Comparison  of  Results  of  Contracted  Wake  Analyses  With 

Goldstein-Lock  Results  and  Experimental  Performance  Results 
for  Model  Rotors  --  0^  =  -8°,  AR  =  13-6,  OR  =  700  fps. 
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Figure  106.  Comparison  cf  Results  of  Contracted  Wake  Analyses  With 

Goldstein-Lock  Results  and  Experimental  Results  for  Model 
Rotors  --  0^  =  -8°,  AR  =  l8.2,  UR  =  525  fps . 
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Figure  1C7.  Concluded. 
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Figure  108.  Sensitivity  of  Predicted  Model  Rotor  ivr f-rsenee 
to  Changes  in  Tip  Vortex  location. 
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Figure  110.  Conparison  of  Results  of  Various  Analyses  With 

Experimental  Performance  Results  for  CH-53A  Rotor. 
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Figure  111.  Sensitivity  of  Predicted  CH-53A  Rotor  Performance 
to  Change  in  Tip  Vortex  Location. 
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MOMENTUM  ANALYSIS 


- GOLDSTE  IN-LOCK  ANALYSIS 

- PRESCRIBED  CLASSICAL  WAKE  ANALYSIS 

- —PRESCRIBED  EXPERIMENTAL  WAKE  ANALYSIS 

(GENERALIZED  WAKE) 


Figure  112.  Comparison  of  CH-53A  Blade  Angle  of  Attack 
Distribution  Predicted  by  Various  Analyses. 
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WAKs  PARAMETERS,  k,  AND  k2 
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Figure  llL.  Variation  of  Cooputed  and  k^  Wake  Parameters 
With  Maximum  Blade  Circulation. 
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APPENDIX  I 

EQUATIONS  FOR  COW  PITTING  THE  VELOCITIES  AND 
DISPLACEMENTS  OF  WAKE  VORTEX  ELEMENTS 


The  velocities  induced  at  an  arbitrary  point  P  by  a  straight  vortex 
element  of  strength  F  and  bounded  by  end  points  A  and  B  (see  Figure  77) 
can  be  computed  using  the  classical  Biot-Savart  Law  (see  Reference  23, 

P.  373)  and  are  given  by  the  following  equations: 

%  “(47)*  [(yp-yA)(*p-zB)  -  (zp-2A)(yp-yB)| 

vyp  1  (4^  ( ^zp  - zaHxp“  xb)  “  ( *p  ~  xa)  (zp  ~  zb)J 

%  1  (47)*  |(xp  _xA)(yp  -yB)-  (*p-yA)  (xp  xb)J 

_  |  (AP  +  BP)/  (AP)  (BP) 

whare  k  '  "R  (AP)(BP)  +•  I  +  J  +  K 

I  -  (*p  -XA)  (Xp  -Xg) 

J  s  (yp-yA){/p  -y*) 

K  8  (zp-za)(2p-z.) 

AP  =  ^/(Xp  -*AfHyP-yAf+  (Zp-ZA)Z  BP  z  ^/\*{>  ~i9?+(zP  ~zB)2 

The  velocities  Vy  ,  and  are  numerically  integrated  using  the 

equations  below  to  determine  the  displacements  of  the  wake  which  occur 
during  a  smell  time  interval,  At 

AXp  2  I  v,pAt 
Ayp  *  I  v>pAt 
AZp  :  l  v,p  Af 


The  summation  signs  in  the  immediately  preceding  equations  denote 
summations  of  the  induced  velocities  induced  at  point  P  by  all  vortex 
elements  in  the  wake.  The  particular  integration  scheme  described  above 
is  based  on  the  assumption  that  the  induced  velocities  remain  essentially 
constant  during  the  time  interval,  At.  Also,  the  length  of  any  vortex 
element  is  allowed  to  vary  as  its  end  points  move. 

Examination  of  the  velocity  equations  given  above  discloses  that  the 
velocity  induced  at  point  P  by  the  straight  vortex  segments  of  which  P 
forms  an  end  point  is  always  zero.  Actually  the  vortex  filament  is  curved 
rather  than  straight.  Some  error  is  thus  introduced  by  the  use  of 
straight  segments,  with  the  size  of  the  error  depending  upon  the  length 
of  the  segment  and  the  curvature  of  the  filament  involved.  To  avoid  this 
error,  the  calculation  of  the  curved  vortex  element  immediately  adjacent 
to  the  point  in  question  was  approximated  by  a  circular  segment,  and  the 
Influence  of  this  circular  segment  was  included  in  the  computational 
program. 


APPENDIX  II 

CONVERGENCE  OF  THE  WAKE  GEOMETRY -BOUND  CIRCULATION  SOLUTION 


As  mentioned  previously,  an  iteratiori  is  required  using  the 
Prescribed  Wake  Program  and  Wake  Geometry  Programs  to  insure  reasonably 
compatible  wake  geometries  and  bound  circulation  distributions.  Because 
of  the  computing  time  required  by  the  Wake  Ceometry  Program,  it  is 
obviously  desirable  to  minimize  the  number  of  passes  through  this  program. 
This  can  be  accomplished  by  (l)  adjusting  the  circulation  distributions 
used  as  input  to  the  Wake  Geometry  Program  so  as  to  anticipate  the  final 
circulation  distribution  answer  as  much  as  possible  (using  past  experience 
as  a  guide)  and  (2)  being  aware  of  the  sensitivity  of  the  final  answer  of 
interest,  namely,  rotor  performance,  to  possible  departures  frcm  the 
ideal,  completely  compatible,  geometry-circulation  solution.  Both 
approaches  were  employed  in  this  study.  The  paragraphs  below  present 
some  results  which  can  be  used  to  evaluate  the  expected  sensitivity  of 
rotor  performance  to  departures  from  the  ideally  converged  solution. 

Rotor  performance  is,  of  course,  given  in  terms  of  integrated  rotor 
thrust  (Ct/ot)  and  torque  (Cq/<t).  Of  these,  rotor  thrust  exhibits  the  most 
sensitivity  to  the  k^  and  kp  wake  geometry  parameters  for  the  tip  vortex. 
Partial  derivatives  relating  C-r/cr  to  k^  and  k2  for  two-bladed  rotors  were 
estimated  by  using  the  Prescribed  Wake  Program  and  varying  k^  and  k2 
independently.  The  results  indicated  the  following  approximate  relation 
for  two-bladed  rotors  : 


ACt/<x  s  ~044  Ak,-0.5Ak2 


(6) 


Now,  if  k^  and  kg  could  be  related  to  the  bound  circulation  distribution, 
one  would  be  able  to  assess  the  probable  impact  of  further  refinements  to 
the  circulation  by  usinc  Equation  (6)  to  compute  an  equivalent  error  in 
Cx/cr  and  comparing  this  with  the  level  of  accuracy  to  which  C7/cris 
desired.  Assuming  that  the  peak  circulation  r on  the  blade  is  the 
characteristic  quantity  determining  the  flew  field  and  hence  k^  and  k2, 
one  can  use  the  computed  wake  geometry  results  that  have  been  obtained 


under  this  contract  for  the  various  operating  conditions  to  assess  the 
rates  of  change  of  k^  and  k2  with  respect  to  r^.  Such  results  are 
shewn  in  Figure  114  for  the  -8-degree  twist  rotor  cases.  Values  are 
plotted  both  as  functions  of  the  specific  values  computed  for  the 

model  rotor  and  a  nondimensional  to  increase  the  utility  of  the 

chart. 

If  C i/cr  is  desired  to  an  accuracy  of  *0.001  (typically  *1^1) »  then 
from  Equation  (4)  this  limits  k^  and  k2  errors  to  *0.0023  and  *0.002, 
respectively.  From  Figure  114  the  corresponding  tolerances  on  r ^  for 
the  model  rotor  are  *19  and  *1.8  ft2/sec,  respectively.  The  large 
tolerance  in  r associated  with  the  error  in  k^  simply  reflects  the 
relative  insensitivity  of  k^  to  (or  thrust)  as  predicted  by  the 

Wake  Geometry  Program.  This  insensitivity  appears  to  be  due  to  the  fact 
that  the  general  dovnwash  induced  by  the  wake  cm  the  tip  vortex  trailed 
by  a  given  blade  tends  to  be  cancelled  by  an  upwash  induced  by  the 
contracted  vortex  trailed  by  the  immediately  preceding  blade. 

All  of  the  geometry-circulation  results  presented  herein  have 
converged  to  within  the  smallest  tolerance  in  I'mfiLX  quoted  above.  Thus, 
the  corresponding  computed  CT/j-  should  be  accurate  to  *0.001. 
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